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ABSTRACT 

With the goal to study the physical and chemical evolution of ices in solar-mass systems, a spectral 
survey is conducted of a sample of 41 low luminosity YSOs (L ~ 0.1 — 10 L Q ) using 5-38 /im Spitzer 
Space Telescope and 3-4 /im ground-based spectra. The sample is complemented with previously 
published Spitzer spectra of background stars and with ISO spectra of well studied massive YSOs 
(L ~ 10 5 L©). This paper focuses on the origin of the prominent absorption features in the 5-8 /im 
spectral region. The long-known 6.0 and 6.85 /im bands are detected toward all sources, with the 
Class 0-type low mass YSOs showing the deepest bands ever observed. In almost all sources the 6.0 
/zm band is deeper, by up to a factor of 3, than expected from the bending mode of pure solid H2O, 
based on the optical depths of the 3.0 /zm stretching and 13 /im libration modes. The depth and shape 
variations of the remaining 5-7 /zm absorption indicate that it consists of 5 independent components, 
which, by comparison to laboratory studies, must be from at least 8 different carriers. Together with 
information from the 3-4 /zm spectra and the additionally detected weak 7.25, 7.40, 9.0, and 9.7 /jm 
features it is argued that overlapping bands of simple species are responsible for much of the absorption 
in the 5-7 /zm region, at abundances of 1-30% for CH 3 OH, 3-8% for NH 3 , 1-5% for HCOOH, -6% for 
H2CO, and ~0.3% for HCOO~ with respect to solid H2O. The 6.85 /zm band likely consists of one or 
two carriers, of which one is less volatile than H2O because its abundance relative to H2O is enhanced 
at lower H2O/T9.7 ratios. It does not survive in the diffuse interstellar medium (ISM), however. The 
similarity of the 6.85 /im bands for YSOs and background stars indicates that its carrier(s) must 
be formed early in the molecular cloud evolution. If an NHj]~ salt is the carrier its abundance with 
respect to solid H2O is typically 7%, and low temperature acid-base chemistry or cosmic ray induced 
reactions must have been involved in its formation. Possible origins are discussed for the carrier of 
an enigmatic, very broad absorption between 5 and 8 /im. It shows large depth variations toward 
both low- and high-mass YSOs. Weak evidence is found that it correlates with temperature tracers. 
Finally, all the phenomena observed for ices toward massive YSOs are also observed toward low mass 
YSOs, indicating that processing of the ices by internal ultraviolet radiation fields is a minor factor 
in the early chemical evolution of the ices. 

Subject headings: infrared: ISM — ISM: molecules — ISM: abundances — stars: formation — infrared: 
stars — astrochemistry 



1 Some of the data presented herein were obtained at the W.M. 
Keck Observatory, which is operated as a scientific partnership 
among the California Institute of Technology, the University of 
California and the National Aeronautics and Space Administra- 
tion. The Observatory was made possible by the generous financial 
support of the W.M. Keck Foundation. 

2 The VLT/ISAAC spectra were obtained at the European 
Southern Observatory, Paranal, Chile, within the observing pro- 
grams 164.1-0605, 69.C-0441, and 272.C-5008 

3 Division of PMA, Mail Code 105-24, California Institute of 
Technology, Pasadena, CA 91125, USA 

4 AURA/NOAO-South, Gemini Science Center, Casilla 603, La 
Serena, Chile 

5 current address: IPAC, NASA Herschel Science Center, Mail 
Code 100-22, California Institute of Technology, Pasadena, CA 
91125, USA (email: aboogert@ipac.caltech.edu) 

6 Division of GPS, Mail Code 150-21, California Institute of 
Technology, Pasadena, CA 91125, USA 

7 Hubble Fellow 

8 Department of Astronomy, University of Maryland, College 
Park, MD 20742, USA 

9 Leiden Observatory, PO Box 9513, 2300 RA Leiden, the 
Netherlands 

10 SRON, PO Box 800, 9700 AV Groningen, the Netherlands 

11 Department of Astronomy, University of Texas at Austin, 1 



1. INTRODUCTION 

The infrared spectra of protostars and obscured 
background stars show prominent absorption features 
at 3.0, 4.25, 4.7, 6.0, 6.85, 9.7, an d 15 urn along with 
a su ite of weaker fea t ures ( e. g. Id'Hendecourt et al.l 
1991 IWhittet et al. I [1996T lGerakines et all |1999|, 



Schutte et ail l999. Gibb et al. 2000, Kca ne et al.ll2001bl 
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Pontoppidan et all l2003al IGibb et aT]|2004 iKnez et al 



20051 IWhittet et al.l 120071 : see iBoogert fe Ehrenfreunc 
20041 for a complete list of features). These are at- 
tributed to absorption in the vibrational modes of 
molecules in ices, except for the 9.7 /xm band which is 
mostly due to silicates. At the low temperatures (T <20 
K) of dense clouds and circum-protostellar environments 
atoms and molecules freeze out rapidly on dust grains. 
Grain surface chemistry efficiently forms new, simple 
speci es, such as H2O and H2CO (e.g. iTielens fe Hagenl 
1982). Complex species (e.g. polyoxymethylene ['POM', 
-(CH2-0) n -], and hexamethylenetetramine ['HMT', 
C6H12N4]) can be formed through the impact of en- 
ergetic photons or cosmic rays on the i ces, as many 
labor at ory studies have sho w n (e.g. ISchutte et al.l 



Gerakincs et al 



1996, 



19931 iBernstein et al.l I1995L iGreenberg et al 



Moore fc HudsonT 



1995 



1998 



Palumbo et al.ll200Ct iMunoz Caro et al.ll2004fh Thus far, 



only the simple species H 2 0, CO, CO2, CH4, CH3OH, 
NH3 and the 13 CO and 13 C02 isotopes have been 
positively identified in the ices toward both low and 
high mass protostars as well as extincted background 
stars whose lines of sight trace quiescent dense cloud 
material. Reasonable evidence exists for the presence 
of HCOOH, OCS and the ions NHj, OCN", and 
HCOO - as well, although their existence is sometimes 
debated because of inaccurate fits with laboratory 
spectra or the lack of multiple bands for independent 
confirmation. The identification of the ions in the 
ices has been partic ularly controversial. The first ion 
claimed was OCN" ([G rim fc Greenbergj|1987l ) and was 
thought to be produced by heavy energetic processing 
of CO:NH 3 ices. Later it was realized that acid-base 
chemistry in a HNCOrNHs ice c ould yield the same 
products (|Novozamskv et al]|2001f) . Acid-base reactions 
are very efficient and occur at temperatures as low 
as 10 K (jRaunier et al.l l2004bf ). The created ions are 
less volatile than neutral species and would be able to 
form an interstellar salt after the other species have 
sublimated. The presence of complex species formed by 
ultraviolet (UV) photon and cosmic ray processing of 
the simple ice s is at least as co n trove rs ial. Many claims 
were made ([Greenberg et al.1 119951: IGibb fc Whittet] 
2002), but the observational evidence is not firm. It is 
crucial to determine the complexity of the ices in the 
circumstellar environment of low mass Young Stellar 
Objects (YSOs), as they may be delivered directly to 
comets and they may ultimately serve as the source of 
volatiles in planets. At sublimation fronts in the warm 
inner regions of disks they are the starting point of a 
complex gas phase chemistry. The relative abundances 
of N-, C-, and O-bearing species coming from the ices 
determine directly the type of chemistry in hot cores 
(or 'h ot corinos' for low mass objects; e.g. ICazaux et al.l 
[2003]) . 

As a step toward better characterizing the molecular 
content of icy grain mantles, this work presents spec- 
tra of a large sample of low mass protostars and back- 
groun d stars obtained w ith the InfraRed Spectrometer 
(IRS; iHouck et al.ll2004ft at the Spitzer Space Telescope 
([Werner et al.l 20041 ) in the 5-35 /im wavelength range. 
These data are complemented with spectra at 2-4 fim 
obtained with ground-based facilities. Such full coverage 
over the near and mid-infrared wavelength ranges is es- 



sential in determining the composition of the ices; most 
species have several vibrational modes and more secure 
identifications can be made when the features are stud- 
ied simultaneously. The identification process is aided by 
studying large samples of sight-lines facilitating correla- 
tion of band strengths and shapes with each other and, if 
known, with physical characteristics along the sight-lines, 
such as luminosity and evolutionary stage of any heating 
source. The great sensitivity of the IRS allows for ob- 
servations of objects with luminosities that are down by 
orders of magnitude with respect to what could be ob- 
served before in the mid-infrared (typically 1 Lq versus 
10 4 L Q ). 

Most data presented in this paper were obtained from 
spectral surveys of nearby molecular clouds and isolated 
dense cores in the context of the Spitzer Legacy Program 
'From Molecular Cores to Planet-Forming Disks' ('c2d'; 
lEvans et al.ll2003T ). Initial results from this program em- 
phasize the importance of thermal processing in the evo- 
lution of the ices, as the ices surrounding the low mass 
YSO HH 46 IRS are mo re processed compared to B5 
IRS1 ([Boogert et al.|[2004l) . Ices towa rd the edge-on disk 
CRBR 2422.8-3423 were presented in iPontoppidan et all 
(2005) and also show signs of heating. Finally, mid- 
infrared spectra of ices toward background s tars tracing 
quiesc ent cloud material were published in IKnez et alJ 
(2005), showing for the first time that the (unknown) 
carrier of the 6.85 /im absorption band is abundant even 
at the coldest conditions away from star formation. 

This paper specifically addresses questions on the 
identification of the 6.0 and 6.85 /im bands. These 
prominent absorption features are still not fully iden- 
tified, despite being detected nearl y 30 years ago with 
the Kuiper Airborne Observatory ([Puetter et al.l fl979T ) 
and commonly observed toward ma ssive YSOs with the 
Infrared Space Ob servatory (ISO; ISchutte et alJ I1996L 
iKeane et al.ll200TrJ ). The peak position of the 6.85 /xm 
band varies dramatically, a nd is thought to be a func- 
tion of the ice temperature (|Keane et al. 2001b). Either 
one carrier with a pronounced temperature dependence 
of its absorption bands, or two independent carriers, one 
much more volatile than the other, must be responsible 
for the 6.85 /xm feature. If the former is the case, the 
amm onium ion (NHj~) is consi dered a promising candi- 
date (|Schutte fc Khannal[2003l ). The 6.0 /im band was 
initially thought to be main ly due to th e bending mode 
of H 2 ([Tielens et al.lll984l : however, see lCox|[l989[ ). but 
analysis of ISO spectra indicated that this is not the 
case toward several massive YSOs ([Schutte et al.l fl996: 
IKeane et ai1l2001tf ). The depth of the 6.0 /jm band is 
in some cases significantly (factor of 2) deeper than that 
expected from the 3.0 /jm H2O stretching mode. Part 
of this 'excess' might be due to a strong vibration of the 
HCOOH molecule. Recent observations of background 
stars show a small excess (<25%), of which half could be 
due to the enhanced b and strength of t he H2O bending 
mode in C02-rich ices ([Knez et al.ll2005[ ). Finally, it was 
claimed that much of the 6.0 /im excess and part of the 
6.85 /im band is due to highly processed ices, i.e. a mix- 
ture of complex species with C-H and O-H bonds pro- 
duced after irradiation of simple ices ([Gibb fc Whittet] 
2002). In this case, one would expect the 6.0 /im excess to 
be enhanced in high radiation environments, such as near 
more evolved protostars. Thus, fundamental questions 
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remain on the identification of the 6.0 and 6.85 /im bands, 
and with the large sample presented in this work the pos- 
sible answers can be further constrain ed. Subsequent pa- 
pers will specifically address the CO2 (jPontoppidan et all 
l2008l) . CH 4 (jOberg et alJl2008h . and NH 3 (S. Bottinelli 
et al., in preparation) species. Further papers in this 
series will investigate the ices toward background stars 
behind large clouds (C. Knez et al., in preparation) and 
isolated cores (A. C. A. Boogert et al., in preparation). 

In |J2]the source sample is presented, while the observa- 
tions and data reduction of the ground-based L-band and 
Spitzer 5-38 /mi spectra are described in SJ3J In if|] the 
continuum level for the multitude of absorption features 
is determined, and subsequently the best value for the 
H2O column density derived, followed by an empirical 
decomposition of the 5-7 /xm absorption complex. The 
positively identified species, or correlations of the compo- 
nents with other observables are discussed in $5j Finally, 
the abundances of the simple species, and constraints on 
the carriers of unidentified components are discussed in 

m 

2. SOURCE SAMPLE 

The source sample is selected based on the presence 
of ice absorption features and consists of a combina- 
tion of known low mass YSOs and new ones identi- 
fied from their Spitzer /IRAC an d MIPS broad-band 
SEDs (Table HI lEvans et al.ll2007f) . For the known ob- 
jects, their mid-infrared (5-20 /im) spectra are presented 
here for the first time, with the e xception of B5 IRS1, 
HH 4 6 IRS (jBoogert et all l2004t iNoriega-Crespo et al.l 
l200l . CRBR 2422.8-3423 (jPontoppidan et alj 2005). 
and L1 489 IRS, IRAS 041 08+2803B, HH 300, and DG 
Tau B (|Watson et al.ll2004fh The YSOs are located in the 
Perseus, Taurus, Serpens, and Corona Australis molec- 
ular cloud complexes, and a number of nearby isolated 
dense cores (Table [T]). YSOs in the Ophiuchus cloud are 
notably absent, because the Spitzer/IRS spectra were not 
available at the time this paper was written. The spec- 
tral energy distributions of the sample of 41 low-mass 
YSOs span a wide range of spectral indices a = —0.25 
to +2.70, with a defined as 

dlog(XF x ) 
dlog(X) 

including all photometry (F\) available between the A = 
2.17 ^m if s -b and from the Two M icron All Sky Sur- 
vey (2MASS; iSkrutskie et "all l200l and the A = 24 
/tin Spitzer/MIPS b and. In the infrare d broad-band 
classification scheme (W ilking et al.| [2001). most objects 
(35 out of 41) fall in the embedded Class 0/1 category 
(a >0.3). The remaini ng 6 objects are F lat-type ob- 
jects (-0.3 < a < 0.3; iGreene et all 11994ft . Although 
most objects fall in the Class 0/1 category, the dispersion 
within it is large. Some objects are envelope-dominated 
Class sources undetected in the 2MASS i£~ s -band (e.g. 
L1455 SMM1, Bl-c), while others are more likely to be 
Flat- type sources with dispersed envelopes and/or face- 
on disks extincted by foreground dust. The latter sources 
were identified by excluding the X s -band flux in the a 
determination, and are marked as such in Table [T] Thus, 
our sample represents a wide range of YSO evolutionary 



stages, and probably source orientation. 

To contrast the possibly processed ices toward YSOs 
with ices toward quiescent regions, published mid- 
infrared spectra of 2 background stars toward the Ser- 
pens and Taurus clouds were included in the analysis 
(jKnez et al.1l2005| ). In addition, 8 massive YSOs and 2 
low mass YSOs observed with the ISO satellite were in- 
cluded in the sample (Table [I]). These serve as a ref- 
erence in the ice feature analysis, because their mid- 
infrared spectra are well-studied at higher spectral res- 
olution than is possible with Spitzer/IRS and because 
the massive YSOs may trace different physical conditions 
that may affect the formation and evolution of the ices. 

3. OBSERVATIONS AND DATA REDUCTION 

Spitzer/IRS spectra were obtained as part of the c2d 
Legacy program (PIDs 172 and 179) as well as a ded- 
icated open time program (PID 20604). In addition, 
several previously p ublished GTO spectra were included 
(| Watson et alj 12004ft . The coordinates and IRS mod- 
ules that were used for each source are listed in Table [TJ 
For all modules the 2-dimensional Basic Calibrated Data 
(BCD) spectral images produced by the SSC pipeline 
version S13.2.0 or later were used as the starting point 
of the reduction. The low resolution 'Short-Low' and 
'Long-Low' modules (SL and LL; R = A/AA = 60 - 120) 
were reduced in a way that is customary for ground- 
based spectra. First, the spectral and spatial dimensions 
were orthogonalized, and then the 2-dimensional images 
of the two nodding positions were subtracted in order 
to remove extended emission. In some cases the back- 
ground emission is highly structured or multiple sources 
are present in the slit, and instead of subtracting nod- 
ding pairs, the background emission was determined from 
neighboring columns in the extraction process. A fixed 
width extraction was performed and the 1-dimcnsional 
spectra were then averaged. Subsequently the spectra 
were divided by spectra of the standard stars HR 2194 
(AO V; PID 1417; AOR keys 0013024512 and 0013024768 
for SL1 and SL2 respectively) and HR 6606 (G9 III; 
PID 1421; AOR key 0013731840 for LL) reduced in the 
same way in order to correct for wavelength-dependent 
slit losses. Spectral feat ures were di v ided o ut using the 
photospheric models of iDecin et alJ (|2004[ L The spec- 
tra of the 'Short-High' and 'Long-High' modules (SH 
and LH; _R~600) were reduced using the c2d pipeline 
(|Lahuis I [2007) . An important difference with the low 
resolution modules is the need for the removal of spectral 
fringes, which is done in the c2d pipeline as well. Many of 
the Spitzer spectr a were complemente d by ground-based 
Keck/NIRSPEC (TMcLean et al.l ll998ft and VLT/ISAAC 
(|Moorwood I I1997D L-band spectra at resolving pow- 
ers of i?=2000 and 600 respectively. These were re- 
duced in a way standard for ground-based long-slit spec- 
tra. Finally, the Spitzer/IRS modules and ground-based 
spectra were multiplied along the flux scale in order 
to match Spitzer/IRAC 3.6, 4.5, 5.8, and 8.0 /im and 
Spitzer/MIPS 24 /an photometry (jEvans et al.ll2007f ). us- 
ing the appropriate filter profiles. 

For the ISO/SWS satellite spectra used in this paper, 
the latest SPD pipeline version 10.1 products were taken 
from the ISO archive and the detector scans were cleaned 
from cosmic ray hits and averaged. The final spectra do 
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TABLE 1 

Source Sample 



Source 


RA a 


Dec a 


Cloud 


Type b 




Obs. ID d 


Module 1 


L-band m 




J2000 


J2000 














L1448 IRS1 


03-25 09 44 


+30-46-21 7 

I 'J U ■ T: U * iJ J- * 1 


Perseus 


low 


0.34 


0005656832 


SL SH LH 

kj -I— J j kJJ- -L j J—J-L J. 


NIRSPEC 


IRAS 03235+3004 


03:26:37.45 


+30:15:27.9 


Perseus 


low 


1.44 


0009835520 


SL,LL 


NIRSPEC 


IRAS 03245+3002 


03:27:39.03 


+30:12:59.3 


Perseus 


low 


2.70 


0006368000 


SL,SH 




L1455 SMM1 


03:27:43.25 


+30:12:28.8 


Perseus 


low 


2.41 


0015917056 


SL,SH,LL1 




RNO 15 


03:27:47.68 


+30:12:04.3 


Perseus 


low 


-0.21 


0005633280 


LL1,SL,SH,LH 


NIRSPEC 


L1455 IRS3 


03:28:00.41 


+30:08:01.2 


Perseus 


low 


0.98 


0015917568 


SL,SH,LL1 




IRAS 03254+3050 


03:28:34.51 


+31:00:51.2 


Perseus 


low 


0.90 


0011827200 


LL1,SL,SH,LH 


NIRSPEC 


IRAS 03271+3013 


03:30:15.16 


+30:23:48.8 


Perseus 


low 


2.06 


0005634304 


LL1,SL,SH,LH 


NIRSPEC 


IRAS 03301+3111 


03:33:12.85 


+31:21:24.2 


Perseus 


low 


0.51 


0005634560 


SL,SH,LH 


NIRSPEC 


Bl-a 


03:33:16.67 


+31:07:55.1 


Perseus 


low 


1.87 


0015918080 


SL,SH,LL1 


NIRSPEC 


Bl-c 


03:33:17.89 


+31:09:31.0 


Perseus 


low 


2.66 


0013460480 


SL,SH,LL1 




Bl-b 


03:33:20.34 


+31:07:21.4 


Perseus 


low 


0.68 


0015916544 


SL,LL 




B5 IRS3 


03:47:05.45 


+32:43:08.5 


Perseus 


low 


0.51 k 


0005635072 


T T -1 (~1 T HTT T TT 

LL1,SL,SH,LH 


NIRSPEC 


B5 IRS1 C 


03:47:41.61 


+32:51:43.8 


Perseus 


low 


0.78 k 


0005635328 


SL,SH,LH 


NIRSPEC 


L1489 IRS f 


04:04:43.07 


+26:18:56.4 


Taurus 


low 


1.10 


0003528960 


SL,SH,LH 


NIRSPEC 


IRAS 04108+2803B f 


04:13:54.72 


+28:11:32.9 


Taurus 


low 


0.90 


0003529472 


SL,SH,LH 


NIRSPEC 


HH 300 f 


04:26:56.30 


+24:43:35.3 


Taurus 


low 


0.79 


0003530752 


SL,SH,LH 


NIRSPEC 


DG Tau B f 


04:27:02.66 


+26:05:30.5 


Taurus 


low 


1.16 


0003540992 


SL,SH,LH 


NIRSPEC 


HH 46 IRS C 


08:25:43.78 


-51:00:35.6 


HH 46 


low 


1.70 


0005638912 


SL,SH,LH 


ISAAC 


IRAS 12553-7651 


12:59:06.63 


-77:07:40.0 


Cha 


low 


0.76 


0009830912 


LL1,SL,SH,LH 




IRAS 13546-3941 


13:57:38.94 


-39:56:00.2 


BHR92 


low 


-0.06 


0005642752 


SL,SH,LH 




IRAS 15398-3359 


15:43:02.26 


-34:09:06.7 


B228 


low 


1.22 


0005828864 


SL,SH,LL1 




Elias 29 g 


16:27:09.42 


-24:37:21.1 


Oph 


low 


0.53 k 


26700814 


SWS01 sp3 




CRBR 2422.8-3423 h 


16:27:24.61 


-24:41:03.3 


Oph 


low 


1.36 


0009346048 


SL,SH,LH 


NIRSPEC 


RNO 91 


16:34:29.32 


-15:47:01.4 


L43 


low 


0.03 


0005650432 


SL,SH,LH 


ISAAC 


IRAS 17081-2721 


17:11:17.28 


-27:25:08.2 


B59 


low 


0.55 


0014893824 


SL,SH,LL1 


NIRSPEC 


SSTc2dJ171122.2-272602 


17:11:22.16 


-27:26:02.3 


B59 


low 


2.26 


0014894336 


SL,LL 




2MASSJ17112317-2724315 


17:11:23.13 


-27:24:32.6 


B59 


low 


2.48 


0014894592 


SL,LL 


NIRSPEC 


EC 74 


18:29:55.72 


+01:14:31.6 


Serpens 


low 


-0.25 


0009407232 


SL,SH,LH 


NIRSPEC 


EC 82 


18:29:56.89 


+01:14:46.5 


Serpens 


low 


0.38 k 


0009407232 


SL,SH,LH 




SVS 4-5 


18:29:57.59 


+01:13:00.6 


Serpens 


low 


1.26 


0009407232 


SL,SH,LH 


ISAAC 


EC 90 


18:29:57.75 


+01:14:05.9 


Serpens 


low 


-0.09 


0009828352 


SL,SH,LH 




EC 92 


18:29:57.88 


+01:12:51.6 


Serpens 


low 


0.91 


0009407232 


SL,SH,LH 


NIRSPEC 


CK4 


18:29:58.21 


+01:15:21.7 


Serpens 


low 


-0.25 


0009407232 


SL,SH,LH 




R CrA IRS 5 


19:01:48.03 


-36:57:21.6 


CrA 


low 


0.98 


0009835264 


SL,SH,LL1 


ISAAC 


HH 100 IRS' 


19:01:50.56 


-36:58:08.9 


CrA 


low 


0.80 


52301106 


SWS01 sp4 




CrA IRS7 A 


19:01:55.32 


-36:57:22.0 


CrA 


low 


2.23 


0009835008 


SL,SH,LH 


ISAAC 


CrA IRS 7 B 


19:01:56.41 


-36:57:28.0 


CrA 


low 


1.63 


0009835008 


SL,SH,LH 


ISAAC 


CrA 1RAS32 


19:02:58.69 


—37:07:34.5 


CrA 


low 


2.15 


0009832192 


nr CTTTT1 

SL,SH,LL1 




L1014 IRS 


21:24:07.51 


+49:59:09.0 


L1014 


low 


1.28 


0012116736 


SL,LL 


NIRSPEC 


IRAS 23238+7401 


23:25:46.65 


+74:17:37.2 


CB 244 


low 


0.95 


0009833728 


SL,SH,LH 


NIRSPEC 


W3 IRS5 1 


02:25:40.54 


+62:05:51.4 




high 


3.53 


42701302 


SWS01 sp3 




MonR2 IRS3 1 


06:07:47.8 


-06:22:55.0 




high 


1.66 


71101712 


SWS01 sp3 




GL989 1 


06:41:10.06 


+09:29:35.8 




high 


0.52 


71602619 


SWS01 sp3 




W33A 1 


18:14:39.44 


-17:52:01.3 




high 


1.92 


32900920 


SWS01 sp4 




GL7009S 1 


18:34:20.91 


-05:59:42.2 




high 


2.52 


15201140 


SWS01 sp3 




GL2136 i 


18:22:26.32 


-13:30:08.2 




high 


1.48 


33000222 


SWS01 sp3 




S140 LRSP 


22:19:18.17 


+63:18:47.6 




high 


1.57 


22002135 


SWS01 sp4 




NGC7538 IRS9 1 


23:14:01.63 


+61:27:20.2 




high 


2.31 


09801532 


SWS01 sp2 




Elias 16 k 


04:39:38.88 


+26:11:26.6 


Taurus 


bg 




0005637632 


SL,SH 




EC 118 k 


18:30:00.62 


+01:15:20.1 


Serpens 


bg 




0011828224 


SL.SH 





a Position used in Spitzcr/IRS observations 

Source type: 'low'— low mass YSO, 'high'— massive YSO, 'bg'— background star 
c Broad-band spectral index as defined in Eg. FTl 

AOR key for Spitzer and TDT number for ISO observations 
e Published previously in Boogert ct al. 2004 
^ Published previously in Watson et al. 2004 
g Published previously in Boogert et al. 2000 
k Published previously in Pontoppidan ct al. 2005 
1 Published previously in Kcanc ct al. 2001 
J Published previously in Knez ct al. 2005 

a. enhanced due to foreground extinction. Exclusion iC s -band flux gives much lower a: —0.16 (Elias 29), —0.02 (B5 IRS1), 0.18 (B5 
IRS3), and 0.38 (EC 82) 

1 Spitzer/IRS modules used: SL=Short-Low (5-14 /mi, R ~ 100), LL— Long-Low (14-34 fim, R ~ 100), SH= Short- High (10-20 /xm, R ~ 600), 
LH^Long-High (20-34 ^m, R ~ 600); ISO SWS modes used (2.3-40 pm): SWS01 speed 1 (R ~ 250), speed 2 (R ~ 250), speed 3 (R ~ 400), 
speed 4 (R ~ 800) 

m Complementary ground-based L-band observations with Keck/NIRSPEC or VLT/ISAAC 
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not show any sig nificant differ e nces w ith respect to the 
data published in lKeane et al.1 (|2001bl ). 

4. RESULTS 

All sources in the sample show numerous absorption 
features on top of rising, flat, or falling 2-35 /im Spec- 
tral Energy Distributions (Fig.Q]). Most of these features 
are due to ices ( N4.2p , and they are often accompanied by 
strong 9.7 and weaker 18 /im silicate absorption bands; 
but toward some sources, the silicate bands are in emis- 
sion (e.g. RNO 15, IRAS 17081-2721). Also, promi- 
nent narrow emission lines of Eb (S[5], S[6], and S[7]) 
are observed toward IRAS 03271+3013, CrA IRAS32, 
and L1455 SMM1. These belong to the most deeply em- 
bedded objects (a > 2; Table [1]) in the sample and the 
line emission is most likely related to mo l ecular outflow 
activ ity (e.g. iNoriega-Crespo et alj 120041 ; iNeufeld et al.1 
120061 ). CrA IRAS32 shows the shock-rel ated 5.340 //m 
emission line of Fe+ ([Neufeld et al.l [2006 ). and a strong 
line at 6.628 /im, perhaps also due to Fe + . These emis- 
sion lines are not further analyzed here, but their pres- 
ence compromises the analysis of the underlying absorp- 
tion features in these particular sources. H I emission 
lines (most notably Br a at 4.052 /im and Pf 7 at 3.740 
/im) are present in many of the ground-based L-band 
spectra, but at the high spectral resolution contamina- 
tion with the ice absorption features is not a problem. In 
sources with strong L-band H 1 emission (e.g. RNO 91) 
the blended Pf a (7.459 /im) and Hu (3 (7.502 /im) lines 
are weakly visible in the IRS spectra (~3% of the con- 
tinuum), affecting the continuum determination for the 
low c ontrast ice features in that region ( §4.2) Obcrg et al. 
120081) . 

4.1. Continuum Determination 

Determining the baseline for the absorption features 
in the spectra of embedded protostars is not trivial. The 
continuum SED underlying the many blended absorp- 
tion features originates from warm dust near the cen- 
tral star and is dependent on poorly characterized pa- 
rameters such as disk and envelo pe mass, size, and in - 
clination, and dust composition (jWhitnev et al.l [2003) . 
The continuum applied to each source and indicated in 
Fig. [1] was derived as follows. Short-ward of 5 /im the 
continuum is much better defined than at longer wave- 
lengths, and therefore the 2-5 and 5-35 /im regions are 
treated separately. For the short wavelength part, con- 
tinuum points selected at 2.3, 3.9, and 5.0 /im are inter- 
polated with a smooth spline function. Inclusion of the 
2.3 /xm portion of the spectrum, or i4T s -band photome- 
try, is essential in choosing the best continuum for the 
3.0 /xm absorption band. The continuum long-ward of 
5 /xm is much less constrained, because the entire 5-32 
/im region is covered by absorption features. Even the 
'emission peak' at 8.0 /im is covered by the long wave- 
length wing of the 6 .0 /im H 2 ice absorption band (e.g. 
iHudgins et al.l fl993h . and the short wavelength wing of 
the 9.7 /xm silicate band. The latter often extends to- 
ward shorter wavelengths compared to the silicate spec- 
tr um observed in the diffuse medium (e.g. toward GCS 
3; lKemper et al.ll2004D . Some sources also show an addi- 
tional depression at 8 /xm due to yet unknown species; the 



best known example of this is the massive YSO W33 A 
dKeane et alJl2001bl : iGibb fc Whittetl [200l . Thus, an 
interpolation over a large wavelength range is required, 
with few constraints. When possible, a low order poly- 
nomial (on a log(XF\) versus log(X) scale) was used to 
interpolate between the 5 and 32 /im continuum points. 
This works best for the most deeply embedded sources 
(e.g. Bl-c, 2MASSJ17112317). For many sources, some 
guidance for the polynomial fit is provided by scaling 
th e silicate spectrum o f the diffuse medium source GCS 
3 ([Kemper et al.|[200l) to the peak optical depth at 9.7 
/xm and by scaling a laboratory H 2 ice spectrum to the 
3.0 /im band. Less embedded sources, i.e. the sources 
with relatively strong near-infrared emission (e.g. RNO 
91, IRAS 03235+3004), show a decreasing continuum 
between 5-10 /im and an emission bump at longer wave- 
lengths. Here, a best effort spline continuum was defined, 
guided by modified black-body curves below 10 /im. For 
th e background stars the optical depth spectra presented 
111 iKnez et al.l {2005) were used. 

4.2. Absorption Features 

The derived continua were used to put the spectra on 
an optical depth scale for the analysis of the wealth of 
absorption features (Fi gs. D] and Et. The w ell known 6.0 
and 6.85 /xm bands fe.g. iKeane et al.l2001bl ) are detected 
in all sources, along with the 9.7 /xm band of silicates 
(either in absorption or emission) and the 3.0 /xm H2O 
and 15 /im CO2 ice bands when available. A number of 
weaker ice features (7.25, 7.40, 7.67, 9.0, and 9.7 /im) are 
detected in several sources as well. Examples in Fig. [3] 
show that these same absorption features are present in 
sources of a very different nature (background stars, mas- 
sive YSOs, Class 0, 1, and flat spectrum low mass YSOs), 
but that the depths and the profiles vary. Very deep 6.0 
and 6.8 /im bands, with peak optical depths of up to 
2, are observed in the Class sources. These are the 
deepest ice bands ever detected in low mass YSOs, com- 
parable to the most obscured massive YSOs (e.g. W33A, 
GL 7009S). 

4.2.1. The H2O Ice Column Density 

H2 O is the most abundant molecule in in terstellar ices 
(e.g. iKeane et all 12001 ri IGibb et al.ll2004D . and strong, 
broad absorption bands are present over much of the ob- 
served spectral range: the O-H stretching mode between 
~2.7-3.6 /im, the O-H bending mode between ~5.4-9 /xm, 
and the libration mode between ~10-30 /xm. Many ab- 
sorption features of other species overlap with these H2O 
bands, especially in the 5-8 /xm region f §4.2.2ll . and in or- 
der to separate them, accurate (10-20%) solid H2O col- 
umn densities are required. The 3.0 /im band provides, 
in principle, the most accurate H2O column density, be- 
cause the O-H stretching mode is very stron g and domi- 
nates over other species ([Tielens et al.|[l984l) . Evidently, 
the ice band must not be saturated, the continuum 
signal-to-noise values must be sufficient, and photome- 
try must be available below the 2.8 /xm atmospheric gap 
for an accurate continuum determination ( §4.1|l . None 
of these requirements are met for the most embedded 
objects in our sample (e.g. Bl-c, 2MASSJ17112317- 
2724315). For sources that do meet the requirements, the 
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Fig. 1. — Observed Spitzcr/IRS and complementary ISO/SWS and ground based L-band spectra and the adopted co ntinua (dotted lines). 
The large dots represent 2MASS and Spitzer/IRAC and MIPS photometry from the c2d catalogs IIEvans et al.1120071) . Sources are ordered 
from top to bottom and left to right in decreasing 9.7 (im absorption band depth. Massive YSOs and background stars are separated in 
the right panel below the white space. 
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Fig. 2. — Optical depth 5-20 fim spectra, ordered from top to bottom and left to right in decreasing 9.7 fim absorption band depth. 
The horizontal dotted line indicates zero optical depth for reference. Massive YSOs and background stars are separated in the right panel 
below the white space. 
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Fig. 3. — Optical depth spectra between wavelengths of 5.2-38 fim (left p anel), and 5.2-8.0 (right panel) of a wide variety of sight-lines . 
From top to bottom: EC 118, a background star behind the Serpens cloud (Kricz ct al. 2005), W 33A, a massive YSO (Kcanc ct al. 2001b), 
and the low mass YSOs Bl-c in the Perseus cloud, L1014 IRS in the isolated core L1014, L1489 IRS in the Taurus cloud, and RNO 91 in 
the L43 core near the Ophiuchus cloud. The low mass YSOs are ordered from top to bottom in decreasing 2-25 /im SED spectral index a 
(Eq. [TJ as indicated in the left panel. Also indicated is the bolometric luminosity. The spectra are scaled to the peak optical depth r max 
in the plotted range indicated near each spectrum. Note that despite the very different nature of the sight-lines, the same ice absorption 
features are observed, although the shapes and depths vary from source to source. 



long wavelength wing that often accompanies the 3.0 /im 
band is not taken into account in the H2O column deter- 
mination. The width of a laboratory pure H2O absorp- 
tion spectrum at the best fitting temperature is assumed 
(Hud ghTs et al.|[T993h . using th e integrated band s trength 
A = 2.0 x 10" 16 cm/molecule (jHagen et al.lll98lD . Thus, 
the calculated column density represents that of absorp- 
tion by pure H2O on small grains. If the wing were 
taken into account, the H2O column density would be 
larger by at most 15% in some sources. The wing, which 
has an integrated optical depth of at most 30% with re- 
spect to the rest of the 3 /im band, is probably caused 
by a combination of scattering on large ice coated grains 
(radius >0.5 um) and absorption by ammonia hy drates 
|Hagen et al.lll983tlDartois fc d'HendecourtJl200l( i. Only 
the former would contribute to the H2O column density, 
but less than expected for small grains because the in- 
tegrated extinction (absorption+scattering) strength per 
molecule is a factor of 2 larger for large grains than it is 
for sm all grains (e.g. Fig. 7 in lDartois fc d'Hendecourtl 
l200l . 

An alternative, and sometimes the only, way to deter- 
mine the H2O column density is by measuring the peak 
strength of the libration mode in the 11 — 14 /im re- 
gion. This band is difficult to isolate in the spectrum 
because it is several /im wide and wedged between the 
strong 9.7 and 18 /im silicate features. However, appli- 
cation of a low order polynomial continuum to HH 46 
IRS (Fig. |4|) and subtract ion of a standard in terstellar 
silicate spectrum (GCS 3; iKemper et al.ll2004|) shows a 
broad feature centered on ^12.3 /im that is attributed to 
the H2O libration mode. Using a laborato ry spectrum of 
a thi n film of pure H 2 ice at T =10 K (Hudgins ct al. 
1993) the depth is in good agreement with the 3.0 /an 



stretching mode, but there is a discrepancy in the peak 
position by about 1 /im (Fig. [5]). The peak position of 
the libration mode is quite sensitive to the shape of the 
grains, a nd for small spherical grai ns (using the formu- 
lation of iBohren fc Huffman~l 1 19831 ) a much better cal- 
culated fit is obtained (Fig. [4]). At such large shift to 
shorter wavelengths, the H2O libration mode overlaps 
with the steep edge of the 9.7 /im silicate band and a 
substructure is introduced re miniscent of the 11.3 /im 
absorption feature reported in lBoogert et al.l (|2004D and 
iKessler-Silacci et al.1 (|200"5h . Finally, this method of de- 
termining the H2O column density works well for deeply 
embedded sources (Fig. [5]), but it fails when a silicate 
emission component is present, filling in the H2O ab- 
sorption (e.g. RNO 91). 

4.2.2. The 5-8 /tm Absorption Complex 

Prominent absorption bands are present at 6.0 and 6.8 
/xm in all sources. In-between these bands (~6.5 /im), 
the optical depth is non-zero and varies relative to the 
6.0 and 6.8 /im peak optical depths for different sources 
(e.g. Fig. [3]). Beyond 7.0 /im, 'featureless' absorption is 
present up to 8 /im that approaches zero optical depth 
at 8 /im in some sources, but continues to be significant 
in others and merges with the 9.7 /im silicate band. Here 
an attempt is made to decompose this complex blend of 
absorption features. Distinct, much weaker bands that 
are present in a subset of sources at 7.25, 7.4, and 7.6 
/im are discussed separately fi j4.2.3p . 

Solid H2O is a major contributor to the optical depth 
in the 5-8 /im region by its O-H bending mode at 6.0 
/im and the overtone of t he libration mode extending be- 
tween 6 and 8 /im (e.g. iDevlin et al.1 1200 If ). For consis- 
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Fig. 4. — Proof of the presence of the 13 /im solid H2O libra- 
tion mode in the source HH 46 IRS. Panel a shows a smooth, 
unbiased polynomial continuum (thick line). The optical depth 
spectrum is compared to the silicate s pectrum of the galac tic cen- 
ter source GCS3 in Panel b (thick line; IKemper et al.l2004D . After 
subtracting GCS3 a prominent absorption feature in found at 10- 
20 /im which is likely due to the strong libration mode of H2O 
(panel c), as is shown by a comparison to a laboratory spectrum 
of pure amorphous H2O ice at T =10 K (thick line; panel d). 
This band is particularly sensitive to grain shape effects, and a 
better fit is obtained with small spherical grains (thick line; panel 
e). The column densities of H2O derived from the 3 and 13 /im 
bands are in good agreement, further validating the correctness of 
the applied method. The ver tical dashed line is the 11.3 (im ab- 
sorption feature identified in IBo ogcrt & Ehrcnfrcund ( 2004|) and 
IKessler-Silacci et al.l d2005t) . 



tency with previous work, although formally incorrect, 
these overlapping bands will be referred to as a single 
mode, the O-H bending mode of H 2 0. Its absorption 
profile depends on the ice structure. In particular, as 
an amorphous ice is heated from T = 1 to 80 K, the 
distinct peak at 6.0 p b ecomes weaker (Hudgin s et al.l 
Il993t iMaldoni et al.lll998h . At higher temperatures, the 
shape remains relatively constant. For almost half of the 
sample the ice temperature could be determined from 
the shape of the 3.0 /jm H 2 stretching mode (Tabled). 
Most other YSOs are deeply embedded, and T — 10 K 
was assum ed. Subsequently a laboratory spectrum of 
pure H 2 (Hudg ins et al.lll993l ) at the determined tem- 
perature was scaled to the column density derived from 
the 3.0 and 13 [im bands f ^4.2.ip and subtracted from 
each optical depth spectrum. The shapes and depths 
of the residual features, and any variations thereof be- 
tween sources, are analyzed. It should be kept in mind, 
however, that some features may still be the result of 
interactions of species with H2O in solid state matrices, 
such as the narrow peak induced in H2 0:C02 mixtures 
( m^l [KneilTal1[20rj5l: lOberg et al.ll2Ml . 

Toward most low mass YSOs, the observed 6.0 /Ltm 
absorption band is deeper than expected from the 3.0 
and 13 /im H2O ice bands (Fig. [6]). Such a discrepancy 
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Fig. 5. — In deeply embedded sources the 13 /im H2O libration 
mode can be used to reliably measure the solid H2O column den- 
sity. 'Silicate free' optical depth spectra were derived in the same 
way as shown in Fig. [4] for the YSOs OA IRS7 B (left panels) 
and Bl-b (right panels). For each source the top panel shows the 
observed spectrum and the polynomial continuum (dotted line), 
and the bottom panel the optical depth spectrum. Note the very 
prominent libration mode and the peculiar absence of silicates to- 
ward Bl-b. Again, small spherical grains of pure H2O at T=10 
K (thick line) fit the peak position of the observed band well, al- 
though excess absorption is visible on the long wavelength side. 
Finally, note that the down-turn above 30 fim is likely real and 
due to the onset of the H2O lattice mode, as the good fits with the 
laboratory data show. 
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Fig. 6. — A selection of 5-8 (im spectra of four low mass YSOs 
on normalized optical depth scales (histogram line). The normal- 
ization factors are indicated for each source. The smooth solid line 
represents a laboratory spectrum of pure, amorphous solid H2O at 
the indicated temperature at the column density derived from the 
3.0 and 13 /im bands. Clearly, in all sources a significant fraction 
of the absorption is not explained by solid H2O, and it varies from 
source-to-source. For the 6.0 (im feature this fraction is indicated 
by the value of Eg , defined in Eq. f2] Note the variations of the 
profiles of the 6.0 and 6.85 fim bands, as well as the presence of 
weak features at 7.25, 7.4, and 7.65 (im. 
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Fig. 7. — Values of Eq (defined in Eq. [2j plotted as a function 
of the solid H2O column density for the entire source sample: low 
mass YSOs (filled circles), massive YSOs (open circles), and back- 
ground stars (stars). For sources near the horizontal dotted line 
the entire 6.0 /im absorption band can be explained by pure H2O 
ice. All other sources have significantly deeper 6.0 fim bands than 
expected from the 3 or 13 fim bands. 

has been reported toward massive YSOs as well (|Coxj 



1989; 



Schu tte et all 119961 : iGibb et al.ll2000t iKeane et all 



2001bft and it is quantified in units of the H 2 column 



density as follows: 



E* = 



Jl562 Tv<1V 



xV(H 2 0) x 1.2 x 10- 17 x 0.60 



(2) 



which is the integrated optical depth of the 6.0 /jm band 
in wavenumber (v) space divided by the H2O column 
density derived from the 3.0 and/or 13 /im bands. The 
factor of 1.2 x 10~ 17 is the integrated band strength of 
the O-H bending mod e in a pure H 2 laboratory ice 
(jGerakines et alj |1995[ ) in units of cm/molecule. The 
latter extends all the way out to 8.5 /im, while the nu- 
merator of Eq. [2] integrates only to 6.4 /xm. Therefore 
the laboratory band strength is reduced by a factor of 
0.60, as derived from the laboratory spectra. Eq varies 
significantly: from >2 in -30% of the YSOs (e.g. RNO 
91, IRAS 03301+3111, EC 92) to 1.0 (i.e. no excess) in 
other YSOs (Tableland Fig. 0). 

The large source-to-source variations of Eq may hold 
clues to the nature of the carrier(s) of the 6.0 /im excess 
absorption. As is illustrated in Fig. these variations 
are used to separate any independent, overlapping ab- 
sorption components in the 5-8 (im range. Initially, a se- 
lection of sources with low values of Eq and sources with 
high values of Eq is averaged, and the resulting spectra 
are scaled and subtracted from each other such that the 
distinct peak at 6.85 /xm is removed. What remains is 
a very broad feature, stretching between 5.8 and 8 /im. 
While the profiles of the 6.0 and 6.85 /xm bands them- 
selves look rather similar for many sources (e.g. like HH 
46 IRS), significant source to source variations are some- 
times observed (Fig.|9|). The peak position of the 6.0 /xm 
component varies between 5.85 and 6.2 /im, and broadens 
as the feature shifts toward longer wavelengths. At the 
same time the 6.85 /xm band shifts from 6.8 to 6.95 /xm. 
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Fig. 8. — Demonstration of the decomposition of the 5-8 
fim absorption complex. In the left panel are shown averaged, 
smoothed, H2 0-subtracted spectra of YSOs with small 6.0 fim 
excess (a), YSOs with strong 6.0 (im excess (b), and the YSO 
IRAS 03301+3111 which shows strongly red-shifted 6.0 and 6.85 
fim bands (c). Applying linear combinations of spectra (a), (b), 
and (c) in the left panel, 5 independent absorption components are 
found and shown in the right panel. Components CI and C3 are 
found by subtracting spectrum (c), multiplied by a factor of 0.63, 
from spectrum (a) [top]. Components C2 and C4 are identical 
to spectrum (c) [middle]. Component C5 is found by subtract- 
ing spectrum (a), multiplied by a factor of 0.69, from spectrum 
(b) [bottom]. These 5 components, in addition to the bending 
mode of H2O, are responsible for the 5—8 (im absorption complex 
in YSOs and background stars. Example fits are shown in Fig. 1101 
Note that C3 and C4 are the previously fo und short- and long - 
wavelength components of the 6.85 fim band (Kcanc et al. 2001b). 



By subtracting extreme spectra from each other, it is 
found that the latter consists of two components centered 
at 6.755 and 6.943 /im with widths of FWHM=0.195 and 
0.292 /im, respectively, simil ar to what was prev iously 
found for massive protostars (|Keane et~alll2001bfl . The 
6.0 /im band consists of two components as well, with 
peak positions of 5.84 and 6.18 /im and FWHM widths 
of 0.26 and 0.40 /xm, respectively. 

It is thus hypothesized that, after H 2 subtraction, 
the 5-8 /im absorption in all sources is composed of 5 
distinct components (Fig. [5]): a set of bands peaking at 
5.84 and 6.18 /xm (CI and C2 hereafter), a set of bands 
peaking at 6.755 and 6.943 /im (C3 and C4), and a broad 
asymmetric feature stretching between 5.8 and 8 /xm and 
peaking near 5.9 /im (C5). Example fits with these com- 
ponents are shown in Fig. [TJJJ The peak optical depths of 
the components in all sources are listed in Tabled The 
indicated error bars are based on statistical uncertainties 
in the spectral data points, except for C5. The accuracy 
of the C5 optical depth is affected by the accuracy of the 
baseline, which at present can only partly be quantified. 
The baseline uncertainty is a combination of the error 
in A(H 2 0) (Table!!]) and the poorly known error in the 
continuum choice ( fl4.ip . For many sources the latter, 
not included in Tabled is likely at most 0.03 on optical 
depth scale. 

4.2.3. Weak Features in the 7-8 fim Wavelength Range 

In addition to solid H 2 and the 5 components dis- 
cussed above, a number of weak features can be discerned 
in the 7-8 /im region of several low mass YSOs, much like 
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TABLE 2 

Fit Parameters (publish electronically) 








Source 


T H2 o a E 6 h t 9 .7 c T 6 .o d r G1 c r C2 f 


T"C3 g 


T C 4 h 


TC5 1 



L1448 IRS1 10 4.35 (1.57) 0.36 (0.04) 

IRAS 03235+3004 10 2.01 (0.31) 1.98 (0.20) 

IRAS 03245+3002 10 1.76 (0.25) 8.11 (0.81) 

L1455 SMM1 10 2.39 (0.37) 5.18 (0.52) 

RNO 15 10* 1.90 (0.19) emission 

L1455 IRS3 10 2.45 (1.00) 0.15 (0.05) 

IRAS 03254+3050 40* 1.75 (0.23) 1.63 (0.16) 

IRAS 03271+3013 10 2.75 (0.63) 1.29 (0.13) 

IRAS 03301+3111 40* 2.61 (0.30) 0.49 (0.05) 

Bl-a 10 1.49 (0.32) 2.35 (0.23) 

Bl-c 10 1.48 (0.28) 7.03 (0.70) 

Bl-b 10 1.30 (0.24) 2.30 (0.23) 

B5 IRS3 10* 1.84 (0.18) 0.32 (0.04) 

B5 IRS1 10* 1.64 (0.21) 1.23 (0.12) 

L1489 IRS 40* 1.77 (0.21) 1.32 (0.13) 

IRAS 04108+2803B 40* 1.61 (0.23) 0.49 (0.05) 

HH 300 10* 2.00 (0.20) 1.09 (0.11) 

DG Tau B 100* 1.87 (0.32) 2.05 (0.20) 

HH 46 IRS 40 2.00 (0.20) 2.23 (0.22) 

IRAS 12553-7651 10 0.99 (0.19) 1.40 (0.14) 

IRAS 13546-3941 10 1.19 (0.12) 0.23 (0.03) 

IRAS 15398-3359 10 2.02 (0.54) 3.32 (0.33) 

Elias 29 40* 1.32 (0.13) 1.44 (0.14) 

CRBR 2422.8-342 40* 1.40 (0.14) 1.68 (0.17) 

RNO 91 40* 1.99 (0.17) 0.80 (0.08) 

IRAS 17081-2721 10* 1.61 (0.16) emission 

SSTc2dJ171122. 2-272602 10 1.58 (0.32) 3.09 (0.31) 

2MASSJ17112317-2724315 10 1.46 (0.24) 4.58 (0.46) 

EC 74 10* 1.75 (0.31) 0.47 (0.05) 

EC 82 10 <0.63 emission 

SVS 4-5 10 2.41 (0.48) 1.61 (0.16) 

EC 90 10 1.01 (0.10) emission 

EC 92 40* 2.28 (0.20) emission 

CK4 10 <1.08 emission 

R CrA IRS 5 40* 1.62 (0.12) 1.02 

HH 100 IRS 10* 2.05 (0.20) 0.69 

CrA IRS7 A 10 0.97 (0.17) 3.02 

CrA IRS7 B 10 1.24 (0.22) 3.44 

CrA IRAS32 10 3.47 (1.24) 2.43 

L1014 IRS 10 1.52 (0.19) 2.08 

IRAS 23238+7401 10 1.17 (0.20) 2.30 

W3 IRS5 10* 1.03 (0.10) 6.28 

MonR2 IRS3 100* 2.41 (0.24) 2.25 

GL989 40* 1.65 (0.08) 0.66 

W33A 40 2.57 (0.64) 6.02 

GL7009S 10 1.98 (0.40) 3.62 

GL2136 100* 1.58 (0.16) 3.02 

S140 IRS1 100* 1.74 (0.17) 1.93 

NGC7538 IRS9 10* 1.58 (0.16) 2.25 

Elias 16 10* 1.48 (0.11) 0.73 

EC 118 10 1.17 (0.12) 1.42 



(0.10) 
(0.07) 
(0.30) 
(0.34) 
(0.25) 
(0.21) 
(0.23) 
(0.63) 
(0.23) 
(0.07) 
(0.60) 
(0.36) 
(0.30) 
(0.19) 
(0.23) 
(0.07) 
(0.14) 



0.10 

1.35 

2.92 

1.92 

0.08 

0.13 

0.32 

0.95 

0.06 

0.73 

2.04 

1.15 

0.12 

0.19 

0.37 

0.22 

0.24 

0.23 

0.75 

0.15 
0.15 
1.41 
0.24 
0.32 
0.39 
0.10 
1.00 
1.30 
0.10 
0.04 
0.68 
0.13 
0.24 
0.08 
0.29 
0.27 
0.54 
0.68 
0.86 
0.62 
0.72 
0.29 
0.20 
0.18 
1.87 
1.17 
0.35 
0.16 
0.48 
0.18 
0.23 



(0.01) 
(0.07) 
(0.17) 
(0.08) 
(0.01) 
(0.01) 
(0.02) 
(0.06) 
(0.01) 
(0.04) 
(0.08) 
(0.04) 
(0.02) 
(0.01) 
(0.02) 
(0.01) 
(0.01) 
(0.02) 
(0.04) 
(0.01) 
(0.01) 
(0.08) 
(0.02) 
(0.03) 
(0.02) 
(0.01) 
(0.05) 
(0.06) 
(0.01) 
(0.01) 
(0.04) 
(0.01) 
(0.02) 
(0.00) 
(0.02) 
(0.03) 
(0.04) 
(0.04) 
(0.04) 
(0.04) 
(0.03) 
(0.01) 
(0.01) 
(0.01) 
(0.03) 
(0.05) 
(0.01) 
(0.01) 
(0.02) 
(0.01) 
(0.02) 



0.00 

0.38 

0.52 

0.55 

0.01 

0.08 

0.14 

0.39 

0.00 

0.33 

0.65 

0.51 

0.01 

0.07 

0.15 

0.10 

0.04 

0.09 

0.34 

0.01 
0.05 
0.54 
0.05 
0.06 
0.13 
0.04 
0.24 
0.35 
0.02 
0.02 
0.31 
0.05 
0.15 
0.03 
0.11 
0.10 
0.08 
0.19 
0.59 
0.36 
0.29 
0.00 
0.01 
0.06 
0.46 
0.38 
0.11 
0.05 
0.12 
0.04 
0.06 



(0.01) 
(0.15) 
(0.40) 
(0.24) 
(0.01) 
(0.04) 
(0.06) 
(0.11) 
(0.00) 
(0.13) 
(0.37) 
(0.21) 
(0.01) 
(0.01) 
(0.06) 
(0.03) 
(0.03) 
(0.02) 
(0.13) 
(0.03) 
(0.03) 
(0.17) 
(0.02) 
(0.05) 
(0.05) 
(0.02) 
(0.14) 
(0.24) 
(0.03) 
(0.02) 
(0.11) 
(0.02) 
(0.04) 
(0.03) 
(0.06) 
(0.05) 
(0.09) 
(0.12) 
(0.29) 
(0.15) 
(0.10) 
(0.04) 
(0.01) 
(0.03) 
(0.20) 
(0.18) 
(0.05) 
(0.02) 
(0.08) 
(0.02) 
(0.05) 



0.00 

0.13 

0.00 

0.17 

0.04 

0.08 

0.04 

0.16 

0.03 

0.12 

0.14 

0.05 

0.04 

0.03 

0.05 

0.06 

0.05 

0.08 

0.17 

0.02 
0.03 
0.15 
0.04 
0.09 
0.07 
0.01 
0.13 
0.04 
0.03 
0.00 
0.09 
0.02 
0.07 
0.04 
0.09 
0.10 
0.03 
0.06 
0.41 
0.21 
0.07 
0.00 
0.02 
0.04 
0.28 
0.01 
0.04 
0.03 
0.07 
0.04 
0.04 



(0.01) 
(0.15) 
(0.40) 
(0.24) 
(0.01) 
(0.04) 
(0.06) 
(0.11) 
(0.00) 
(0.13) 
(0.37) 
(0.21) 
(0.01) 
(0.01) 
(0.06) 
(0.03) 
(0.03) 
(0.02) 
(0.13) 
(0.03) 
(0.03) 
(0.17) 
(0.02) 
(0.05) 
(0.05) 
(0.02) 
(0.14) 
(0.24) 
(0.03) 
(0.02) 
(0.11) 
(0.02) 
(0.04) 
(0.03) 
(0.06) 
(0.05) 
(0.09) 
(0.12) 
(0.29) 
(0.15) 
(0.10) 
(0.04) 
(0.01) 
(0.03) 
(0.20) 
(0.18) 
(0.05) 
(0.02) 
(0.08) 
(0.02) 
(0.05) 



0.02 

0.40 

1.19 

0.86 

0.03 

0.04 

0.18 

0.33 

0.00 

0.39 

1.01 

0.52 

0.06 

0.07 

0.19 

0.10 

0.08 

0.08 

0.36 

0.09 
0.05 
0.55 
0.05 
0.13 
0.13 
0.07 
0.48 
0.74 
0.05 
0.02 
0.28 
0.06 
0.10 
0.01 
0.13 
0.11 
0.26 
0.35 
0.34 
0.34 
0.38 
0.14 
0.06 
0.10 
0.74 
0.72 
0.19 
0.05 
0.27 
0.07 
0.20 



(0.01) 
(0.07) 
(0.16) 
(0.16) 
(0.01) 
(0.03) 
(0.02) 
(0.13) 
(0.01) 
(0.04) 
(0.06) 
(0.04) 
(0.02) 
(0.02) 
(0.02) 
(0.01) 
(0.02) 
(0.02) 
(0.02) 
(0.03) 
(0.01) 
(0.03) 
(0.02) 
(0.02) 
(0.01) 
(0.01) 
(0.04) 
(0.04) 
(0.02) 
(0.01) 
(0.02) 
(0.01) 
(0.02) 
(0.02) 
(0.01) 
(0.03) 
(0.02) 
(0.02) 
(0.22) 
(0.04) 
(0.05) 
(0.01) 
(0.02) 
(0.01) 
(0.05) 
(0.04) 
(0.02) 
(0.01) 
(0.01) 
(0.01) 
(0.03) 



0.03 

0.13 

0.55 

0.25 

0.03 

0.02 

0.09 

0.07 

0.03 

0.29 

0.64 

0.35 

0.05 

0.03 

0.08 

0.06 

0.06 

0.10 

0.18 

0.07 
0.02 
0.31 
0.02 
0.11 
0.06 
0.03 
0.30 
0.42 
0.02 
0.01 
0.21 
0.03 
0.02 
0.01 
0.08 
0.09 
0.35 
0.31 
0.31 
0.16 
0.16 
0.16 
0.22 
0.07 
0.51 
0.54 
0.14 
0.10 
0.15 
0.03 
0.10 



(0.01) 
(0.07) 
(0.16) 
(0.16) 
(0.01) 
(0.03) 
(0.02) 
(0.13) 
(0.01) 
(0.04) 
(0.06) 
(0.04) 
(0.02) 
(0.02) 
(0.02) 
(0.01) 
(0.02) 
(0.02) 
(0.02) 
(0.03) 
(0.01) 
(0.03) 
(0.02) 
(0.02) 
(0.01) 
(0.01) 
(0.04) 
(0.04) 
(0.02) 
(0.01) 
(0.02) 
(0.01) 
(0.02) 
(0.02) 
(0.01) 
(0.03) 
(0.02) 
(0.02) 
(0.22) 
(0.04) 
(0.05) 
(0.01) 
(0.02) 
(0.01) 
(0.05) 
(0.04) 
(0.02) 
(0.01) 
(0.01) 
(0.01) 
(0.03) 



0.08 

0.39 

0.84 

0.83 

0.02 

0.00 

0.04 

0.26 

0.02 

0.03 

0.24 

0.00 

0.02 

0.01 

0.06 

0.01 

0.04 

0.01 

0.10 

0.00 
0.00 
0.40 
0.00 
0.00 
0.10 
0.01 
0.19 
0.37 
0.01 
0.00 
0.17 
0.00 
0.02 
0.00 
0.00 
0.00 
0.00 
0.04 
0.00 
0.00 
0.00 
0.00 
0.09 
0.01 
0.67 
0.37 
0.06 
0.00 
0.04 
0.00 
0.00 



(0.01) 
(0.06) 
(0.13) 
(0.07) 
(0.01) 
(0.03) 
(0.01) 
(0.04) 
(0.01) 
(0.05) 
(0.13) 
(0.07) 
(0.02) 
(0.01) 
(0.02) 
(0.02) 
(0.01) 
(0.01) 
(0.02) 
(0.02) 
(0.01) 
(0.09) 
(0.02) 
(0.02) 
(0.02) 
(0.01) 
(0.06) 
(0.07) 
(0.01) 
(0.01) 
(0.03) 
(0.01) 
(0.01) 
(0.02) 
(0.01) 
(0.02) 
(0.04) 
(0.04) 
(0.05) 
(0.02) 
(0.05) 
(0.02) 
(0.01) 
(0.01) 
(0.07) 
(0.06) 
(0.01) 
(0.01) 
(0.02) 
(0.02) 
(0.01) 



Note. — Uncertainties in parentheses based on statistical errors in the spectra only, unless noted otherwise below. 
a Temperature of pure H2O laboratory ice (^udgins ct al. 1993) assumed. Values indicated with an asterisk are most accurate because of the availability 
of good quality L-band spectra. 

6.0 fiui excess as denned in Eq.[2] Uncertainty includes uncertainty in N(H20). 
c peak optical depth at 9.7 /iin, without correction for underlying emission. Uncertainty includes 10% of T9.7 to account for errors in the continuum 
determination. 

peak optical depth at 6.0 /im, including H2O absorption. 
e peak optical depth component CI 

peak optical depth component C2 
g peak optical depth component C3 

peak optical depth component C4 

1 peak optical depth component C5. Uncertainty includes uncertainty in A^l^O) as this affects the baseline level. The continuum uncertainty is not 
included. It is likely on the order of 0.03 for most sources. 
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Fig. 9. — Residual absorption in the 5-8 fim region for a selection 
of sources, after subtraction of pure H2O and component C5. The 
vertical dashed lines at 6.0 and 6.85 /an are plotted to facilitate 
comparing the positions of the features in the different sources. 
Note that the plotted sources represent the most extreme cases of 
our sample; most sources resemble the spectrum of HH 46 IRS. The 
spectra are normalized on optical depth scale, and sorted from top 
to bottom in increasing wavelength of the 6 fim residual. Note that 
at the same time the 6.85 fim feature shifts to longer wavelengths 
as well. 
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Fig. 10. — Examples of fits to the 5-8 fim absorption complex 
with the five components derived in Fig. [8] as well as H2O. In the 
left and right panels two different YSOs are shown. The dots are 
the observed data points, and the thin line through them is the 
sum of all components above. 
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Fig. 11. — Optical depth 7-8 (im Spitzer/IRS spectra of sources 
showing the 7.25, 7.40, and 7.67 (im features. The top source, 
W33A, is a massive YSO of which a higher resolution (ij=800) 
ISO/SWS spectrum is shown. The features are labeled with their 
(possible) identifications. 



those identified toward massive protostars ([Schutte et al.l 
1999). These features lie on top of a curved 'continuum' 
created by the wings of the prominent 6 and 6.85 /im 
features. A high order polynomial (up to 6 th ) had to 
be used to fit a local baseline to typically these wave- 
length regions: 7.00-7.14, 7.3-7.35, 7.47-7.50, and 7.80- 
7.95 /im. Examples of the resulting features are shown 
in Fig. llli where they are compared to the higher resolu- 
tion (R = 800) ISO/SWS spectrum of the massive YSO 
W33A. T he 7.25 [im feature is likely due to the HCOOH 
molecule (< j5.1.2p . T he feature at 7.40 /xm is possibly due 
to HCOO- C 35X51 IScinTtte et a.l.lfl999h. while th e 7.67 
/jm feature is identified with CH4 (jBoogert et al.lll996l : 
lOberg et al.ll200l . 

4.2.4. Substructures in the 9.1 \iva Silicates Absorption 

Band 



Several sources show substructures at 9.0 and 9.7 /im 
within the strong silicate absorption band (Fig. [T2J), Pre- 
vious work indicated that these can be attributed to the 
umbrella mode of solid NH3 a nd the C-0 stretching mod e 
of CH3OH (e.g. lLacv et al.lfl998t iSkinner et al7fl992]) . 
To extract these and any other ice absorption features in 
the 8-12 /zm range, the silicate absorption and possible 
blended emission needs to be accurately modeled. This is 
deferred to a future publication (S. Bottinclli, in prep.), 
and here a local 4 th order polynomial continuum is fit- 
ted to the wavelength regions 8.25-8.75, 9.23-9.37, and 
9.98-10.4 /xm to extract the features (Fig. [TSJ) . Column 
densities of CH3OH and NH3 derived from these and 
other features are further discussed in § §5.1.11 and !5.1.4l 



CONSTRAINING THE ORIGIN OF THE 5-8 
ABSORPTION FEATURES 



Numerous carriers have been proposed to account for 
the various feature s in the 5-8 urn wa velength region 
(|Tielens et al.l 11984 iKeane et al.lboOlbl ). Here new ob- 
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9.0 9.5 
A (/im) 

Fig. 12. — Features of solid NH3 and CH3OH in the spectra of 
the low mass YSOs HH 46 IRS (left panels) and IRAS 15398-3359 
(right). In the top panels for each source the dashed line indicates 
a local polynomial continuum, which was used to derive the optical 
depth plots in the lower panels. 

servational constraints to possible carriers are discussed, 
using the present large source sample and wavelength 
coverage. First, the contributions of the most securely 
identified (CH3OH, NH3) or otherwise well constrained 
species (H 2 CO, HCOOH, HCOCT) to the observed Cl- 
C5 components are discussed ( §5.ip . Subsequently, the 
origin of the remaining absorption is constrained by cor- 
relating it with a number of other observables ( £|5.2|) . 

5.1. Absorption by Known Species 



5.1.1. CH3OH 

Solid CH3OH was previously unambiguously de- 
tected by its C— H stretching mode at 3.54 /im (e.g. 
iGrim et al]|1991h and i ts C— O stretching mode at 9.7 
/im (Skin ner et aHfl992h . In our sample, a Gaussian was 
fitted to the 9.7 /im feature on an optical depth scale and 
a CH3OH detection is claimed if the integrated optical 
depth has a larger than 3c significance and if the peak po- 
sition and FWHM are similar to previous detections, i.e. 
within 9.68-9.76 and 0.25-0.36 /an respectively. The 3.54 
/zm band, when available, was treated similarly. Thus, 
CH3OH was positively identified in 12 low mass YSOs 
(Table [3j Fig. [13]). Of those, 6 are based on the 9.7 
/im band only. Column densities were calculated using 
integrated band strengths of 5.6xl0 -18 and 1.6xl0~ 17 
cm/molecule for the 3.54 and 9.7 fim bands, respectively 
(iKerkhof et al.lll999T). I n agreement with previous stud- 
ies (|Dartois et all 119991: iPontoppidan et al J 120039) , the 
CH3OH column varies considerably with respect to solid 
H2O: from low upper limits of a few percent (3a) in sev- 
eral sources, to detections of up to 30% in a few others. 
In the 5-8 /im region, the C— H deformation mode of 
CH3OH overlaps with the C3 component, and, at the 
derived column densities, contributes typically 5-20% to 
its strength, but up to 40% for SVS 4-5 and GL7009S. 




5.1.2. HCOOH 

The main vibrational modes of formic acid (HCOOH) 
occur in the 3-4 /im region (O-H and C-H stretch), at 5.85 



Fig. 13.— Abundances of solid CH3OH (panel a), HCOOH 
(b), and NH^~ (c) in percentage of H2O as a function of the 9.7 /im 
silicate band optical depth. The solid CH3OH abundance is clearly 
more sight-line dependent than the HCOOH abundance. Sources 
with the largest relative CH3OH column densities are labeled in 
the left panel. Open symbols with arrows refer to 3<r upper limits. 
Error bars are la. The NH^~ abundance in the right panel was 
calculated assuming both the C3 and C4 components are due to 
NH^~, but after subtracting the contributions from CH3OH and 
H 2 CO. 



/im (C=0 stretch), at 7.25 /im (C-H deformation), and at 
8.1 /im (C-0 stretch). Previously, the 7.25 /im band was 
detected toward massive YSOs, while the 5.85 /im band, 
in additi on to H2O, could ex plain most of the 6 /im ab- 
sorption (jSchutte et all 1 19991 ). The 7.25 /im band offers 
the most secure identification, even though its intrinsic 
strength (A = 2.6 xlO -18 cm/molecule) is a factor of ^25 
weaker than the 5.85 /im band. It is least blended with 
other features and it is unambiguously detected (> 3c) 
toward 9 low mass YSOs (Fig. [Til Table 131 W ith the 
aid of laboratory analogs (|Bisschop et al. Il2007t ). the ob- 
served strength of the 7.25 /im feature is used to deter- 
mine the contribution of HCOOH to the 6.0 /im band. 
Pure HCOOH, and mixtures with H2O show 7.25 /im 
features that are at least a factor of 2 wider than the 
observed band (FWHM ^15 cm" 1 ), and that peak at 
too sh ort wavelengths . ISchutte et al.l (|1999f ) and re- 
cently iBisschop et al. I (|2007T > note that the band nar- 
rows and shifts to longer wavelengths in tertiary mix- 
tures, specifically in H 2 O:CH 3 OH:HCOOH=100:40:12 
(Fig. [T4"|) . Such a mixture would be consistent with the 
detections and upper limits of solid CH3OH (typically 
7V(CH 3 OH)> iV(HCOOH) toward YSOs; Table©. For 
the purpose of this analysis the laboratory and interstel- 
lar data are matched by scaling their integrated 7.25 /im 
optical depths. For these sources, the HCOOH column 
densities are 1-5% with respect to solid H2O (Fig. [T3|) and 
50-100% of component CI can be explained by the C=0 
stretching mode of HCOOH (Fig. [TJJ. This wide range 
encompasses the large uncertainty in the band strength: 
the laboratory 5.85 /im band is stronger by 60% in H2O 
mixtures compared to pure HCOOH. For many sources 
in which the 7.25 /im band is not detected (< 3cr), a 
stronger constraint on the HCOOH column density is 
provided by the strength of the CI feature. Relative 
HCOOH abundances of a few % are derived by using the 
C=0 band strength of the H20-rich mixtures. Finally, in 
most cases neither the series of shallow, broad HCOOH 
bands in the 3-4 /im O-H and C-H stretching mode re- 
gion, nor a band at 8.1 /im, blended with the edge of 
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TABLE 3 
Column Densities of the Ices 



Source 


AT(H 2 0) 


Ar(HCOOH) c 


A r (CH 3 OH) 


^(NHJ) 6 




10 18 cm" 2 


% H 2 


% H a O 


% H 2 


L1448 IRSl 


47 (0 161 


<12.1 


<14.9 


13.9 (4.6) 


TRAS 03235+3004 


1 4 48 (2 26l b 


<2.7 


4.2 (1.2) 


2.2 (0.9) 


IKAb Uo24o+oU02 


oy.ol (o.ooj 


<1.2 


<9.8 


5.5 (0.7) 


L1455 bMMl 


18.21 (2.82)" 


3.3 (0.8) 


<13.5 


7.1 (1.6) 


RNO 15 


0.69 (0.06) 


<5.8 


<5.0 d 


13.5 (2.4) 


L1455 IRS3 


0.92 (0.37) b 


<15.1 


<12.5 


8.9 (5.5) 


IRAS 03254+3050 


3.66 (0.47) 


<3.2 


<4.6 


8.2 (1.1) 


IRAS 03271+3013 


7.69 (1.76) b 


<2.5 


<5.6 


5.8 (3.3) 


IRAS 03301+3111 


0.40 (0.04) 


<4.9 


<15.1 


8.1 (3.2) 


Bl-a 


10.39 (2.26) b 


<2.8 


<1.9 


7.7 (0.7) 


Bl-c 


29.55 (5.65) b 


1.2 (0.3) 


<7.1 


6.7 (0.3) 


Bl-b 


17.67 (3.20) b 


3.1 (0.5) 


11.2 (0.7) 


2.9 (0.4) 


B5 IRS3 


1.01 (0.09) 


<4.4 


<8.1 


15.3 (3.6) 


B5 IRSl a 


2.26 (0.28) 


<2.6 


<3.7 


5.2 (1.3) 


L1489 IRS 


4.26 (0.51) 


<2.9 


4.9 (1.5) d 


6.2 (0.8) 


IRAS 04108+2803B 


2.87 (0.40) 


<3.1 


<3.5 


6.9 (0.9) 


xlrl OUU 


z.oy yv.zo) 


^9 ^ 




o.o (i.i ) 


no Tan R 


2 29 (0 391 


<3 .4 


<5.7 


9.3 (1.8) 


1414 A(\ TR^ a 


7 7Q (f\ 77\^ 


9 7 fn 71 


O.O ^U.O ) 


p. Q /n 41 


li\/\o lzoOo-(DOl 


o ofi /n ^f^b 
z.ys (^u.oDj 






HA ^ 1 ^ 


T T5 A Q 1 Qt^/lfl QO/11 
Irx/Ao 1004D-oy41 








q r (r\ 7^ 
O.O ^U. / ) 


1HAS 15398-3359 


14.79 (3.95) 


1.9 (0.3) 


10. o (0.8j 


3.9 (0.3) 


Elias 29 a 


3.04 (0.30) 


<11.3 


<4.9 d 


3.2 (1.0) 


CHBH 2422.8-o42 a 


4.19 (0.41) 


<1.0 


<9.3 


6.7 (0.7) 


jxiNvj yi 


4. ZD (^U.oOJ 






0.4 ^U.OJ 




i qi f n 1 q "\ 
l.ol l^U.lo^ 






y.y (i-oj 


QQTnOfl T1 71 1 OO O OVOfiflO 


io.y4 yZ.oA ) 


o.U (U.DJ 




« A /r\ ^\ 
0.4t ^U.OJ 


01\A A QQ 11 71 1 0'?1 7 070/I Q1 £ 


iy.4y yo.zv) 


z.o (u.oj 




« 7 { C\ A \ 
0./ (^ / 


T7Y" 1 7/1 
I -/ L ( 4 


1 n7 /'n ie\ 
l.lw ^U.lo) 


9 


<.y .o 


O 7 {1 Q\ 
y. / \^-o) 


EC 82 a 


SO (() 07) 


<2.5 


<14.2 


10.2 (6.3) 


SYS 4.5a 


5 65 (1 1 31 


4.6 (1.0) 


25 2 ("^ 5^ 

— ■■ .' ■ — yij.iJJ 


3.6 (0.5) 


EC 90 a 


1.69 (0.16) 


<4.1 


6.8 (1.6) 


4.4 (0.9) 


EC 92 a 


1.69 (0.14) 


<5.3 


11.7 (3.5) d 


2.9 (1.8) 


CK4 


<1.56 








R CrA IRS 5 


3.58 (0.26) 


<4.1 


6.6 (1.6) d 


5.1 (0.4) 


HH 100 IRS a 


2.45 (0.24) 


<2.4 


<9.7 d 


11.6 (2.1) 


CrA IRS7 A 


10.89 (1.92) b 


<0.6 


<3.8 


7.0 (0.4) 


CrA IRS 7 B 


11.01 (1.97) b 


<1.5 


6.8 (0.3) 


5.8 (0.2) 


CrA IRAS32 


5.26 (1.88) b 


<9.8 


<18.1 


15.9 (7.9) 


J_iXUJL^t 1JAO 


7 ifi (o QI \b 


^ A (A ^1 
U.I ^l.OJ 


O.J- ^u.o^ 


70 fi nl 


TRAS 2S2384-7401 


1 2 95 (2 26l b 


<1.5 


<3.6 


4.8 (0.7) 


W3 IRS5 a 


5.65 (0.56) 


<0.4 


<2.8 d 


7.0 (0.4) 


MonR2 IRS3 a 


1.59 (0.15) 


<2.7 


<5.6 d 


25.9 (2.2) 


GL989 a 


2.24 (0.10) 


<2.2 


3.1 (0.1) 


8.4 (0.7) 


W33A a 


12.57 (3.14) 


5.2 (0.7) 


14.7 (3.6) d 


8.1 (0.8) 


GL7009S a 


11.31 (2.26) 


2.5 (0.7) 


31.3 (7.8) d 


4.9 (0.7) 


GL2136 a 


4.57 (0.45) 


5.1 (1.2) 


8.5 (2.1) d 


7.2 (0.9) 


S140 IRSl a 


1.95 (0.19) 


<2.6 


<3.0 d 


9.6 (1.0) 


NGC7538 IRS9 a 


6.41 (0.64) 


<2.0 


7.5 (1.8) d 


5.7 (0.4) 


Elias 16 a 


2.26 (0.16) 


<1.4 


<2.3 


5.2 (0.5) 


EC 118 a 


3.57 (0.35) 


<1.4 


<4.7 


11.6 (1.4) 



Note. — Uncertainties (1 cr) arc indicated in brackets and upper limits arc of 3 er signifi- 
cance. 

a Value obtained from or comparable to previous work or references therein: iDartois et all 
fl999l (CH 3 OH: W3 3A, GL2136, W3 IRS 5, S140 IRSl, NGC 7538 IRS9, MonR2 IRS3, HH 
100 IRS, GL7009S'). lBoo£ert et al.ll2000bl fCH.^OH: Eli as 291. iBoogert et al.|[200l fH^O and 
CH3OH: B5 IR Sl and HH 46 IRS) , iBrooke et al.lfl999l (H 2 0: EC 90V lEiroa fc Hodapplfl989l 
(H 2 0: EC 82), |Keane et al.|[2001bi (H 2 0: W33A, GL2 136, W3 IRS5, S 140 IRSl, NGC7538 
IRS9, MonR2 IRS3, H H 100 IRS, GL989, GL70 09S), iKnez et al.|[2005l (H 2 and CH 3 OH: 
Elias 16 and EC 118), iPontoppidan et aL1f200^ (H 2 and CH 3 OH in EC 92, SVS 4-5), 
IPontoppidan et al.|[2005l (H 2 and CH 3 OH: CRBR 2422.8-3423). 

13 fim H 2 libration mode used for 7V(H 2 0) determination (3 /.im band in other cases). 
c Detections and 3-cr upper limits of iV(HCOOH) based on 7.25 ^m C-H deformation mode. 
Sources for which the 5.85 fim C— O stretching mode provides tighter constraints arc indicated 
with '<'. 

The 3.53 ^m band provides a better or comparable constraint to 7V(CH30H) than does the 
9.7 f-iin band. 

c Assuming both C3 and C4 components arc due to NHJ, which is still a matter of debate 
( fl6.1fr ; integrated band strength A — 4.4 X 10 -17 cm/molcculc assumed (Schuttc Sz Khannal 
2003). 
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Fig. 14. — Comparison of laboratory solid HCOOH and interstel- 
lar spectra. The left panels show the region of the C=0 stretch- 
ing mode, the right panels the C-H deformation mode. From top 
to bottom four YSOs are shown with confirmed 7.25 ^m absorp- 
tion band detections (indicated with vertical dotted lines in the 
right panels). The histograms in the left panels represent com- 
ponent CI of the interstellar spectra, i.e. after subtracting solid 
H2O and components C2 and C5 defined in i|4,2,2l The dashed 
lines represent the spectra before C2 and C5 subtraction. In each 
panel the thick smooth line represents the laboratory ice mixture 
HCOOH:H2O:CH 3 O H=12:100:40 at T = 1 5 K after removal of the 
H2O bending mode {Bisschop ct al. 2007). The laboratory spec- 
tra are scaled such that the integrated optical depth of the 7.25 fim 
band matches the observations. Note that the laboratory spectrum 
of the 7.25 fim band is wider than shown in Bisschop et al. (2007). 
This is due to a different baseline choise at the steep edge of the 
strong C-H deformation mode of CH3OH in the laboratory data. 



the prominent 9.7 /im silicate feature, provide stronger 
constraints to the HCOOH column density. 



5.1.3. H 2 CO 



As w as noted for massive YSOs (e.g. iKeane et al.l 
l2001bl ). the CI component shows a sharper edge on the 
short wavelength side than can be explained by HCOOH 
alone. The C=0 stretching mode of H2CO is a viable 
candidate, and, assuming an inte grated band strengt h 
of A = 9.6 x 10" 18 cm /molecule (jSchutte et all 11993( 1. 
substantially better fits to component CI are obtained 
by adding typically an abundance of 6% of H2CO with 
respect to H2O (Fig. [T5|) . Accurate H2CO abundances 
and uncertainty estimates are hard to establish due to the 
strong blend with HCOOH. This abundance estimate is 
consistent with the non-detection of the C-H stretching 
bands at 3.34, 3.47, and 3.54 fim and leads to a contri- 
bution of 10-35% from H2CO to component CI. At the 
quoted typical abundance, the C-H bending mode vibra- 
tion of H2CO at 6.68 fim contributes <15% to compo- 
nent C3 (assuming A — 3.9 x 10~ 18 cm/molecule from 
ISchutte et aTlll993f) . 



required and accurate column densities will be presented 
elsewhere (S. Bottinelli et al., in prep.). In this work, 
NH3 column densities accurate to within a factor of a 
few are sufficient. The strength of the 9.0 /jm umbrella 
mode is used to calculate the strength of the 6.16 /jm 
N-H bending mode and its contribution to component 
C2. The inflection seen at 9.0 fim in several sources is 
fitted with a Gaussian on optical depth scale and a de- 
tection of NH3 is claimed if the measured peak position 
and FWHM width are within 8.91-9.05 fim and 0.20-0.45 
fim, respectively. Thus, NH3 is positively identified in 17 
low mass YSOs. Column densities were calculated using 
an integrated band strength of 1.3 x 10~ 17 cm/m olecule 
for NH 3 mixed in an H 2 ice (|Kerkhof et al.lll999f ), lead- 
ing to values of typically 3-8% with respect to H2O. At 
these abundances, 10-50% of the C2 component is due 
to solid NH 3 . 

5.1.5. HCOO- 

The formate ion (HCOO~) is an interesting candidate 
to consider because its C=0 stretching mode at 6.33 fim 
overlaps with the C2 component, its C-H deformation 
mode coincides with the 7.40 band detected in several 
sight-lines (Fig. [TTJ, and because it is che mically related 
to th e HCOOH species discussed above (jSchutte et al.l 
1999). Its identification is tentative, however, since ac- 
etaldehyde (CH3HCO) has a feature at 7.40 fim as well, 
but not near the C2 component. If HCOO - were the 
carrier of the 7.40 fim band, its column density would be 
^0.3% with respect to solid H 2 0, an order of magnitude 
less than HCOOH, and its contribution to component C2 
would be < 20%. 

5.1.6. Composite Laboratory Spectrum 

To illustrate that the independently derived column 
densities of solid HCOOH, HCOO~, H 2 CO, NH 3 , and 
CH3OH arc indeed consistent with the 5-7 fim observa- 
tions, laboratory spectra of these ices are scaled and sub- 
sequently added (Fig. [15]). While no attempt is made to 
match the band profiles in detail (profiles are highly de- 
pendent on interactions in the ice matrix and a detailed 
analysis awaits further laboratory work), the agreement 
is encouraging. A considerable fraction of the interstellar 
5-8 fim absorption consists of a complex blend of features 
from these simple ices. In particular, they explain most 
of the observed strength of component CI (50-100%), 
much of C2 (30-70%) and some of C3 (<30%). At the 
same time, Fig. [15] also shows that component C5 and 
most of the prominent 6.85 fim absorption band (C3 and 
C4) are not accounted for. NH|" is a candidate for the lat- 
ter ( §6.1|) . but its column density cannot be determined 
independently and thus this species is not included in 
this plot. 

5.2. Correlations 



5.1.4. NH 3 

The abundance of solid NH3 has long been a matter 
of debate, because all of its absorption bands overlap 
with much stronger bands of H2O and silicates. A com- 
bined in-depth analysis of all these features is therefore 



To obtain further insight into the origin of the Cl- 
C5 components, their strengths were normalized to the 
H2O column density for the entire sample and plotted in 
Fig. [16] As a measure of the degree of correlation the ra- 
tio of the mean values and standard deviations were cal- 
culated (excluding upper limits), yielding values of 1.97, 
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6.5 7.0 
A (/tim) 

Fig. 15. — Simple species contributing to the 6 and 6.85 fim 
bands of the low mass YSO HH 46 IRS, after H2O subtraction 



(iV(H 2 O)=82.8xl0 1 



Laboratory spectra of (from top to 
bottom) solid CH3OH, CH 4 , NH 3 , HCOO", H 2 CO, and HCOOH 
were scaled to column densities derived independently elsewhere 
in the spectrum: 5.6% (CH3OH), 2.2% (CH 4 ), 12.7% (NH 3 ), 0.4% 
(HCOO"), 6.0% (H 2 CO), and 2.0% (HCOOH) with respect to 
H2O. These spectra were subsequently added to the underlying, 
unidentified C5 absorption component (dashed line) and the result 
is shown as a solid thick line. Clearly, the prominent 6.85 fim band 
is not explained by these particular species. 



1.70, 3.61, 1.35, and 0.87 for C1-C5 respectively. Thus, of 
all components, C3 correlates best with H2O, followed by 
CI. The C2 and C4 components show more scatter, and, 
most interestingly, their strengths are largest at the low- 
est H2O abundances. Component C5 correlates poorest 
with H2O, especially when taking into account the tight 
3cr upper limits at the bottom end of the distribution. 
While, of all components, C5 is most susceptible to errors 
in the continuum determination, the C5/iV(H20) varia- 
tion is too large to be explained by that alone f £|4.2.2|) . 
Therefore, the large variations in the 6.0 /im excess (Eq; 
Fig. [7]) can be ascribed primarily to variations in compo- 
nent C5, and much less to CI. 

The different degrees of correlation of the C1-C5 com- 
ponents with H2O indicate that they have different ori- 
gins. To further elucidate their nature, their strengths 
were correlated with a number of tracers of the chemical 
or physical processing history of the ices. A well-known 
tracer is the rati o of broad and na rrow 4.7 //m CO ice 
components (e.g. IChiar et al.lll998l) . often referred to as 
'polar' and 'apolar' ices due to the mixing of CO with 
high and low dipole moment species, respectively. The 
apolar ices have a much lower sublimation temperature 
than the polar ices (T su b — 20 versus 90 K) and thus their 
abundance ratio traces thermal processing. Higher tem- 
peratures are traced by the 15 fim CO2 bending mode. 
Characteristic double peaks arc visible in the crystalline 
ice s produced after heatin g C02:H 2 mixtures to T > 50 
K (|Gerakines et alJll999f ). and the fraction of crystalline 
CO2 ice so derived for many source s in our sample is 
presented in lPon toppidan ct al. (2008). Actual dust tem- 
peratures can be derived from far-infrared colors, but to 
trace the coldest dust components (T <20 K) photom- 
etry at wavelengths > 100 /im is required. For massive 
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Fig. 16. — Peak optical depths of the 5-8 /im components (Fig. 
[8}, normalized to A r (H20), plotted as a function of iV(H20). Error 
bars are 1<t, but open symbols with arrows attached refer to 3cr 
upper limits. At the top right of each panel the ratio of the mean 
over the standard deviation of the sample (excluding upper limits) 
is indicated. Component C3, i.e. the short-wavelength component 
of the 6.85 fim feature, correlates significantly better with iV(H20) 
than the other components do. Of the components that produce 
the 6.0 /im excess, CI correlates better with H2O than does C5. 
Note the increased strengths of C2 and C4 at low iV(H20). 



YSOs this was provided by ISO/LWS spectra, and the 
far-infrared color was f ound to correlate with indicators 
of thermal processing (|Boogert et al.ll2000aD . Unfortu- 
nately, such photometry is not available for most of our 
sample of low mass YSOs. Tracers of ultraviolet radi- 
ation and cosmic ray fluxes are even harder to obtain. 
Traditionally, the strength of the 4.62 /im absorption 
band has been used for this, since the carrier, an 'XCN'- 
bearing species, was thought to be pro duced by the ef- 
fects of UV photons or cosmic rays (e.g. iPendleton et all 
1999). Recent work shows, however, that the potential 
carrier OCN - can be produced by heating an HNCO- 
containing ice formed at cold conditions by grain surface 
chemistry (jvan Broekhuizen et al.ll2004T ). Thus, in this 
case the strength of the 4.62 /im band is in fact another 
tracer of thermal evolution. 

The prominent 6.85 /jm absorption band consists of 
independent short and long-wavelength components (C3 
and C4; Fig. [8|), whose integrated optical depth ratio 
varies considerably (C4/C3=0.5-4). This essentially re- 
flects abundance variations of the carrier of the C4 com- 
ponent, because C4 correlates more poorly with the H2O 
column density than does the C3 component (Fig. IT6|) . 
The C4/C3 ratio is plotted against the observables dis- 
cussed above, and while all of the observables show a 
large spread at C4/C3<1.2, sources with stronger C4 
components all have low solid H2O and CO abundances 
(panels (a) and (b) of Fig. [TT|). These include the low 
mass YSOs IRAS 03301+3111 and DG Tau B, and the 
massive YSOs Mon R2 IRS3 and S140 IRS1. The weak- 
ness of the ice bands in these particular sources hampers 
the correlation study with some of the above-mentioned 
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t. . C4/t. . C3 
int ' lint 

Fig. 17. — The ratio of the integrated strengths of components 
C4 and C3 (red and blue components of the 6.85 /im band) plotted 
against possible ice or dust temperature tracers. Along the y— axis 
are plotted: (panel a) a measure of the solid H2O abundance, 
(panel b) a measure of the solid CO abundance, (panel c) the 
polar/apolar solid CO ratio, (panel d) the fraction of crystalline 
CO2 measured from the 15 /im bending mode, and (panel e) the 
relative abundance of the carrier of the 4.62 /im feature, OCN - . 
Low mass YSOs are indicated by filled dots, massive YSOs by 
circles, and background stars by star symbols. Panels (a) and (b) 
show that sources with prominent C4 components have low ice 
(H2O and CO) abundances. 

observables: the ratio of the polar/apolar solid CO com- 
ponents and the fraction of crystalline CO2 cannot be 
determined if these bands are not or weakly detected. 
Similarly, only a handful of OCN - abundance determina- 
tions are available for these sources. It is thus concluded 
that the C4/C3 ratio is enhanced in environments that 
have low ice abundances. T his is in agreement w ith ear- 
lier work on massive YSOs (|Keane "etaLll2001bl ). which 
showed an increase of the ratio in lines of sight with 
warmer dust temperatures. Potential carriers of the C4 
and C3 components are discussed in £16.11 

The strength of the enigmatic C5 component varies 
widely for both low and high mass YSOs. Toward mas- 
sive YSOs, a trend of increasing strength of the C5 com- 
ponent with the OCN~ abunda nce (4.62 /im feature) was 
claimed (|Gibb fc W hittet 2002), but in the current sam- 
ple this trend is rather weak and is limited by the small 
number of sources with OCN - detections (panel (e) in 
Fig. [T5J) , A stronger trend is observed with the ratio 
of the polar and apolar CO components: sources with 
strong C5 components show larger polar/apolar CO ra- 
tios (panel (c) in Fig. [T5|) . No correlations are observed 
with the H2O, CO, and crystalline CO2 abundances. Po- 
tential carriers are further discussed in 

6. DISCUSSION 



Long-standing key questions in astrochemistry concern 
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the composition of the interstellar ices and the degree to 
which the composition of the ices is altered by heating 
and by the impact of energetic photons and cosmic rays. 
The 3-30 /im spectra presented in this paper extend stud- 
ies addressing these questions from the previously well 
analyzed massive YSOs to low mass YSOs, as well as to 
background stars tracing quiescent cloud material. The 
5-8 /xm region exhibits a very rich absorption spectrum, 
which, apart from the O-H bending mode of H2O, con- 
sists of at least 5 independent components. While most 
of component CI, much of C2 and part of component C3 
is due to the simple species CH3OH, HCOOH, NH3, and 
H2CO, the origin of a substantial amount of absorption 
is still unclear, and this is discussed in more detail below. 

6.1. Constraints on the Carrier (s) of Components C3 
and C4 (the 6.85 /im Band) 

The 6.85 /im band was decomposed into 'blue' and 'red' 
components (C3 and C4). While C3 correlates tightly 
with the H2O abundance, the relative strength of the 
C4 component increases at the lowest H2O abundances 
(Fig. [Trj)) . This could be explained by the formation of 
a new species as the ices are heated and sublimate or it 
could be a pre-existing species that is less volatile than 
H2O whose absorption band (C4) prevails at lower H2O 
abundances. In the first scenario heating must play a 
role, and in the second it may. The detection of the 
C4 component in all YSOs and toward background stars 
favors the second scenario, in which the carrier is not 
formed by heating for its carrier exists already in the 
earliest stages of chemical evolution, i.e. in quiescent 
clouds. Higher dust temperatures in the vicinity of YSOs 
lead to sublimation and lower H2O and CO abundances 
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(Fig. [T7|) and the abundance of the carrier of the C3 
component is reduced more rapidly than that of the C4 
component. 

Does the carrier of the C4 component survive in the 
diffuse interstellar medium? A band at 6.85 /im is de- 
tected in lines of sight toward the Gal actic Center (GC) , 
tracing a large column of diffuse dust (|Chiar et al.l l2000) . 
However, its depth relative to the 9.7 /im silicate band 
(t"6.85/ t '9.7 ~ 0.014) is a factor of 5 less compared to the 
C4 component in the YSOs (factor of 10 compared to 
C3). The weak 6.85 /im band observed toward the GC 
sources is attributed to the C-H as ymmetric deforma - 
tion mode of aliphatic hydrocarbons (|Chiar et al.l f2000). 
Such aliphatic hydrocarbons have strong C-H stretching 
modes with distinct peaks at 3.38, 3.42, and 3.48 /im. 
A band at 3.47 /im is commonly detected toward YSOs 
(|Brooke et al.| [l999), but it does not show these distinct 
peaks and its depth relative to the 6.85 /im band is a 
factor of 3 weaker compared to the GC sources. Thus, 
although the carrier of the C4 component observed in our 
sample is (probably) less volatile than ices, it does not 
survive in the diffuse ISM. Conversely, the hydrocarbons 
observed in the diffuse ISM are not the carrier of the C4 
component toward YSOs. 

In a scenario proposed by ISchutte fc Khannal (|2003f ) 
both the C3 and C4 components arise from an NH^- 
containing ice. Its abundance would be typically 7% 
with respect to H 2 (Tableland Fig. [13|), assuming 
the integrated band strength of 4.4 x 10 -17 cm/molecule 
measured in ISchutte fc Khannal (|2003h . NHj ions can 
be produced either by acid-base reactio ns at low tem- 
peratu re (e.g. from NH 3 :HNCO at 10 K; iRaunier et al.l 
l2004bf l or by ultraviolet or cosmic ray irradiation of sim- 
ple mixtures (e.g. H 2 0:C0 2 :NH 3 :02). The NH| band 
that is produced shifts to longer wavelengths upon heat- 
ing, i.e. from the peak wavelength of the C3 component 
to the peak wavelength of the C4 component. After ad- 
ditional energetic processing and H 2 sublimation a salt 
is formed. Salt s do not sublimate until temperatures 
of -200 K (e.g. ISchutte fc Khannal 120031 ). significantly 
higher than the H 2 sublimation temperature of 90 K 
in space. The processed laboratory spectra show absorp- 
tion at 6.0 /im due to anions (NO^, NO^T, and HCOg"), 
some of which overlaps with the C2 components centered 
at ~6.15 /im that are observed in the YSOs with strong 
C4 components (Fig. [9j §6.3| . While these properties 
are consistent with the observations of YSOs, the pres- 
ence of the C4 component toward all YSOs and back- 
ground stars is hard to e xplain. In the experiments of 
ISchutte fc Khannal (|2003h an NH^ temperature of al- 
most 200 K is required to match the peak position of this 
component. Such a high temperature can be excluded 
for many sources in the sample. Perhaps the impact of 
cosmic rays is responsible for the shift in peak position. 
Subsequent heating by nearby protostars leads to the ob- 
served increase in the C4/C3 ratio and eventually to the 
formation of a salt. 

6.2. Potential Carrier (s) of Component C5 

The enigmatic C5 component that stretches from 5.8 to 
~8.0 /xm is not due to NH 3 , H 2 CO, HCOOH, HCOO~ , or 
CH3OH, because these species do contribute to compo- 



nents C1-C4 fi j5.ip . Also, no correlation is found between 
the strengths of the C5 and the 6.85 /im red (C4) compo- 
nents. Therefore, these have unrelated origins. Heating 
might be involved in the formation of the carrier(s) of the 
C5 component (Fig. [181 £j5T2j) . Previously, three different 
origins were considered in the literature. 

First, it should be noted that the shape of the average 
C5 component resembles that of the spectrum of solid 
H 2 (Fig. [T0| . It is particularly reminiscent of heated 
H 2 (T > 120 K), because of the weakness of the 6.0 
/im peak with respe ct to the long wavelength wing (e.g. 
iMaldoni et al.lll998f ). While the shape of the 3 /im band 
indicates that heated ices are present along several lines 
of sight (e.g., RNO 91), its depth is, by the very defini- 
tion of the C5 component, insufficient. Possibly, due to 
the effects of light scattering at the shorter wavelengths, 
a more direct path to the warm inner regions of the YSO 
is observed at 6 /im compared to 3 /im. Indeed, Monte 
Carlo radiative transfer models of a circumstcllar disk 
show that at inclination angles of 70 — 74° the depth of 
the 6.0 /im H 2 band is e nhanced by a factor of 10 with 
respect to that at 3.0 /im ()Pontoppidan et aT]|2005l ). At 
higher inclinations, the enhancement is a factor of 3-4, 
while at lower inclinations the effect is negligible. Thus, 
the C5 component may be attributed to heated H 2 ice 
for some sources in our sample; for e xample, RNO 91 
is known to possess an inclined disk (jWeintraub et al.l 
|1994|) and has one of the most prominent C5 compo- 
nents. At the other extreme, the weakest C5 compo- 
nents are observed toward background stars (Table [2]). 
Indeed, in these lines of sight the effect of scattering on 
the relative H 2 band depths is expected to be small due 
to the uniformity of the foreground clouds. For many 
YSOs with C5 absorption, however, independent obser- 
vational evidence indicates that the scattering scenario 
is unlikely. For example, the CH3OH column density 
derived from the 3.53 /im C-H stretching mode agrees 
well with tha t derived from the 9.7 /im C-H bending 
mode (iU ISchutte et all 11991 . Also, for high mass 
YSOs the 4.25 and 15 /im bands provided solid C0 2 col- 
umn d ensities that were in ex cellent agreement with each 
other (|Gerakines et al.| [T999) , and in one lin e of sight also 
with t he 2.7 /im C0 2 combination mode (jKeane et alJ 
l2001aj ). Thus, for these YSOs other explanations are re- 
quired. 

A second explanation of the C5 component 
could be the ble n d of negative ions suggested in 
ISchutte fc Khannal (|2003l ). Ammonium (NH|) and 
the anions HCO^, NOg", and NO^ are produced 
after UV photolysis or cosmic ray bombardment of 
H 2 0/C0 2 /NH 3 /0 2 ice mixtures. NH^ has a strong 
band at 6.85 /im (i.e. components C3 and C4 in this 
work), but the bands of the anions overlap such that 
only one very broad feature is observed in the 5-8 /im 
range. While one may expect that in this scenario the 
C3/C4 and C5 components correlate well with each 
other, the observations show they do not. Notably, the 
background stars have strong 6.85 /im components and 
little C5 absorption. The NH^ ion can be produced 
after warming up a simple ice mixture too, however, 
while the HCO^" , NO^, NQ7 anions ma y need more 
energetic input (jSchutte fc Kha nna 2003). This may 
explain the large variation in the strength of the C5 
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component. 

Third, extensive processing of simple ices by UV light 
or cosmic rays produc es an organic refractory residue 
(|Greenberg et al.l [1995) , whose blended O-H and C-H 
stretching and bending modes produce a featur e simi- 
lar to the C5 component (|Gibb fe Whittetll2002f ). Both 
this scenario and the blend of ions discussed above re- 
quire a correlation of the strength of the C5 component 
with independent indicators of the presence of strong 
UV or cosmic ray fluxes, which are currently unavail- 
able. Also, both the organic residue and the salts are 
less volatile than the ices — the y sublimate at tempera- 
tures <500 K in the laboratory JSchutte fc Khannal l200l 



iGreenberg fc d'Hendecourtl 119851 : T.„>,=180 K for H 2 
ice). Therefore, one may expect C5 to be the main ab- 
sorption component in the 5-8 /xm region toward some 
YSOs where all H2O ices have evaporated. Although in 
several sources the C5 component is the dominant con- 
tributor, it always is observed together with H2O (except 
possibly toward L1448 IRS1). Some H2O may be trapped 
in the salt or organic residue at high temperatures, and 
some may be present in unprocessed foreground mate- 
rial. It is also important to note that, similar to C4 
( q6.ip . the C5 component is not observed in the diffuse 
ISM, confirming that their carrier(s) are not as refractory 
as silicates (T su b ^1500 K). 

Finally it is noted that the C5 component may have 
a different origin in different lines of sight. An in-depth 
analysis of the physical conditions in specific lines of sight 
is required to discriminate between the effects of disk 
orientation and the various ice processing mechanisms. 

6.3. Potential Carrier (s) of Component C2 

The C2 component spans the 6.0-6.4 /im wavelength 
range and is likely due to at least three different carriers. 
The N-H bending mode of NH3 at 6.15 /im contributes 
10-50% to this band ( §5.1.4j) . Another likely contribu- 
tor is H2O. Although the spectrum of pure, amorphous 
H2O was subtracted from the 5-8 /im region before deter- 
mining the various component strengths and profiles, the 
spectrum of interstellar H2O may deviate from that of a 
pure, amorphous H2O ice. The O-H bending mode of 
H2O shifts to longer wavelengths, becomes significantly 
narrower, and its peak strength increases if the hydro- 
gen bonding network in bulk amorph ous H2O ice is bro- 
ken by the presence of other species (Ehrcnfr eund et al.l 
119961 : lOberg et al.ll2007h . C0 2 is a realistic species, be- 
cause it is abundant (15-40% with respect to H 2 0) and 
a significant fraction is known to be mixed w ith H2O 
(|Pontoppidan et al.ll20MlGerakines et al.lll999f >. At suf- 
ficiently high CO2 concentrations, the H2O is located in 
small clusters, leading to monomer, dimer, and multi- 
mer bondings that produce a strong peak at 6.12 /Ltm. 
The optical depth of the peak has been calculated for 
the YSO B5 IRS1 (lOberg et all [20071) and the back- 
ground star Elias 16 (jKnez et al.l 120051 ) assuming the 
H2O ice along these lines of sight consists of two com- 
ponents: pure H2O and the mixture H20:C02=2:1. It 
was found that when 25% of the H2O is present in the 
latter, the 3 and 6 fim bands of H2O and the 15 /im 
band of CO2 are fitted simultaneously. Comparison with 
the C2 component in these sources shows that 20-50% 



could be due to this effect. Thus the short wavelength 
side of the C2 component is both due to H20:C02 mix- 
tures and due to NH3. For the long- wavelength side, the 
C=0 stretching mode of the HCOO - anion was previ- 
ously considered, contributing <20% to the C2 strength 
( §5.1.5|) . Other anions, such as NO2", NO^, and HCO^ 
(i j6.ip . may contribute as well. These particular anions 
may be produced in the ener getic processing of simple 
ices (jSchutte fc Khannall2003T ). along with NHj. Indeed, 
one of the most interesting observational aspects of the 
C2 component is that it becomes more pronounced for 
sources with stronger C4 components and lower H2O 
abundances (Figs. l9l and [T6|). Overall the two compo- 
nents do not correlate very well, however, which is consis- 
tent with the result that three carriers of different nature 
contribute to the C2 component. 

6.4. Astrochemical Implications 

CH3OH shows larger abundance variations with re- 
spect to H2O than other solid state species discussed 
in this work (HCOOH, H 2 CO, NH 3 ). Tight 3cr upper 
limits of a few % are reported for a number of YSOs, 
but detections above the 10% level for others (Table [3] 
and Fig. [T3|). The solid HCOOH abundance is typically 
a factor of 2-3 lower than CH3OH and does not show 
such large variations. If both CH3OH and HCOOH are 
formed from CO by grain surface reactions, the lower 
HCOOH abundances are understood by the low atomic 
O/H ratio in the accreted gas, and by the competition 
for O to react to CO2 rather t han with the formyl rad - 
ical (HCO) to form HCOOH (jTielens fc Hagen lTi~982T ). 
Indeed, large solid CO2 abundances (15-40% with re- 
spect to H 2 Q) are observ ed for a wide variety of YSOs 
(jPontoppidan et aI1l2008h . Grain surface chemistry does 
indeed appear a more likely formation mechanism than 
UV processing of the ices. First, while processing of sim- 
ple ices by UV photons leads to CH3OH formation, this 
process is not efficient enough to e xplain the enhance d 
abundances toward massive YSOs dDartois et al. 1999), 
while for low mass YSOs the UV flux is likely too low 
regardless. Second, HCOOH and CH3OH are detected 
even in the most obscured sight-lines in the sample. In 
fact, enhanced CH3OH abundances are detected exclu- 
sively in some of the most embedded YSOs: the low 
mass YSOs SVS 4-5, Bl-b, and IRA S 15398-3359 and th e 
massive YSOs W33A and GL7009S (|Dartois et alJl999h . 
SVS 4-5 is known to tr ace the envelope of the fo reground 
Class object SMM4 (jPontoppidan et al.ll200l . On the 
other hand, toward a few background stars tracing qui- 
escent cloud material a low upper li mit to the CH3OH 
abundance was determined fTable 131 IChiar et al.l ll996: 
iKnez et aT1l2005| ). Thus it is suggested that, in the grain 
surface chemistry scenario, the CH3OH abundance is en- 
hanced at a time during the collapse (of some clouds) 
when the atomic O needed to form H2O and CO2 is ex- 
hausted, while still sufficient CO and H are available to 
form CH3OH ijOberg et al.ll20M) . 

While grain surface chemistry likely forms the simple 
neutral species, different mechanisms are required to pro- 
duce the ions possibly responsible for several of the com- 
ponents in the 5-8 /im range. If the observed C3 and 
C4 components are indeed due to NH4", its abundance 
is typically 7% with respect to H 2 (Fig. [13]). The C4 
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component is observed even in the coldest sight-lines. In- 
deed, low temp erature acid-base chem istry can efficiently 
produce NH4 ()Raunier et alj l2004bl ). The abundance 
of negative counter ions needs to be addressed. While 
the C5 component could h ide the necessary negative ions 
(jSchutte fc Khannall2003h . its strength does not correlate 
with the C3/C4 components. Further investigations on 
the presence of ions in the ices are crucial to understand 
the ice chemistry. Their physical properties are differ- 
ent from neutral species. Ions have higher sublimation 
temperatures, and thus the salts that are formed after 
the volatile species (e.g. H 2 and CO2) have sublimated 
are able to survive at lower extinctions, and then are 
subjected to harsher conditions, forming more complex 
species. For example, the biologically interesting species 
urea (H2NCONH2) is formed after the UV irradiation of 
a NH|OCN" salt (jRaunier et al.ll2004al) . 

7. CONCLUSIONS AND FUTURE WORK 

The present work extends the study of ices over the full 
3-20 /im wavelength range from the previously well stud- 
ied massive YSOs (~ 10 5 L Q ) to low mass YSOs (~1L Q ). 
The following new insights were obtained: 

• Absorption features are detected at 6.0 and 6.85 
/im in all sources, including high- and low-mass 
YSOs and background stars tracing quiescent cloud 
material. An empirical decomposition shows that 
the 5-7 /im absorption complex consists, in addition 
to the H2O bending mode, of a combination of at 
least 5 independent absorption components. 

• In a subset of sources additional weak features are 
detected at 7.25. 7.40, and at 9.0 and 9.7 fxm on 
top of the prominent Si-0 stretching band of sili- 
cates. These features are associated with CH3OH, 
HCOOH, NH 3 , and possibly HCOO" and indicate 
abundances of 1-30%, 1-5%, 3-8%, and 0.3% with 
respect to H 2 0, respectively. The large source- 
to-source solid CH3OH abundance variations are 
likely a result of the conditions at the time of grain 
surface chemistry. 

• Component CI (5.7-6.0 /im) is mostly explained 
by solid HCOOH and H2CO at abundances with 
respect to solid H 2 of typically 1-5% and ~6%, 
respectively. 

• Component C2 (6.0-6.4 /im) also likely arises from 
a blend of several species. Solid NH3 can account 
for 10-50% of the absorption. A similar amount can 
be attributed to absorption by monomers, dimers, 
and small multimers of H2O mixed with a substan- 
tial amount of CO2. The long- wavelength side of 
C2 is potentially due to anions produced by acid- 
base chemistry (e.g. HCOO - ) or energetic process- 
ing. Indeed, a weak correlation with C4 suggests 
the carriers of these bands (possibly salts) are re- 
lated. 

• Components C3 and C4, the 6.85 /im band, show 
the same characteristics for low mass YSOs and 
background st ars as was previously found for mas- 
sive YSOs bv iKeane etall (|2001bh . Their ratio is 



empirically found to be dependent on dust temper- 
ature. The carrier of C4 is likely less volatile than 
that of C3 and than H2O, but it is not observed in 
the diffuse ISM. These characteristics are also con- 
sistent with both components being due to NH4" 
(ammonium ion). The detection of strong 6.85 /im 
bands toward deeply embedded YSOs and back- 
ground stars requires a production at low temper- 
ature (acid-base chemistry or processing by cosmic 
rays) , given the lack of heating sources and stellar 
UV fields. 

• The origin of the very broad component C5 (5-8 
/im) is least understood. It is quite strong in several 
low and high mass YSOs, absent in others and so 
far undetected toward background stars. It is pos- 
sibly related to thermal processing, as a weak cor- 
relation with the ratio of the polar/apolar CO com- 
ponents is observed and its shape resembles that of 
warm H2 O. A blend of absorption b y ions, as pro- 
posed by (jSchutte &: Khannal [2003) . which would 
require both heating and energetic processing of 
the ices, cannot be excluded. The latter could also 
lead to an organic residue, whose shape resembles 
that of the C5 component as well. 

• Weak correlations are found between the absorp- 
tion components in the 5-8 /im range and tracers of 
physical conditions. A more thorough understand- 
ing of the conditions, the source geometry and local 
radiation fields in specific lines of sight is required 
to further constrain the nature of the carriers (in 
particular for components C2-C5) and the impor- 
tance of thermal and energetic processing. To what 
degree does scattering play a role? What are the 
dust temperatures, cosmic ray fluxes, UV radia- 
tion fields and time scales as a function of distance 
along the line of sight? Measurements of source 
fluxes above wavelengths of 100 /im with the Her- 
schel Space Observatory will be valuable, as well as 
ground-based measurements of the 4.62 /im band 
of OCN~ in more lines of sight. In addition, fur- 
ther laboratory work, in particular on the effect of 
cosmic rays on the ices, especially for the identifica- 
tion of components C3 and C4, is required. Finally, 
study of the 5-8 /im features toward a larger sam- 
ple of background stars is crucial to make further 
progress in the identification and processing history 
of interstellar ices. 
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